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Abstract 

We study individual eigenstates of quantized area-preserving maps on the 2-torus which are 
classically chaotic. In order to analyze their semiclassical behavior, we use the Bargmann-Husimi 
representations for quantum states, as well as their stellar parametrization, which encodes states 
through a minimal set of points in phase space (the constellation of zeros of the Husimi density) . 
We rigorously prove that a semiclassical uniform distribution of Husimi densities on the torus 
entails a similar equidistribution for the corresponding constellations. We deduce from this 
property a universal behavior for the phase patterns of chaotic Bargmann eigenfunctions, which 
is reminiscent of the WKB approximation for eigenstates of integrable systems (though in a 
weaker sense) . In order to obtain more precise information on "chaotic eigenconstellations" , we 
then model their properties by ensembles of random states, generalizing former results on the 2- 
sphere to the torus geometry. This approach yields statistical predictions for the constellations, 
which fit quite well the chaotic data. We finally observe that specific dynamical information, e.g. 
the presence of high peaks (like scars) in Husimi densities, can be recovered from the knowledge of 
a few long-wavelength Fourier coefficients, which therefore appear as valuable order parameters 
at the level of individual chaotic eigenfunctions. 



1 Introduction 

We present a detailed exploration of the eigenfunctions of certain quantum maps, corre- 
sponding to classical chaotic maps defined over a torus phase space. Those eigenfunctions 
will be described in phase space by means of the Bargmann-Husimi-stellar family of 
representations. This study is intended to contribute towards understanding the general 
issue of the semiclassical behavior of individual eigenfunctions, which is a basic and largely 
unsolved problem when the classical dynamics is chaotic. 

When this dynamics is integrable, we may consider that the semiclassical behavior of 
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the eigenfunctions is known: it follows a WKB Ansatz 

V^-^^^e'^^/^, ^^0, (1) 

3 

where Sj are the (finitely many) branches of the classical action at energy E and Aj are 
certain invariant 1/2-densities. Thus in principle (disregarding fine points and technical- 
ities), the eigenfunctions and all derived quantities arc computable in the semiclassical 
regime (e.g., the discrete eigenvalues are yielded by EBK quantization formulae). 

By contrast, no definite behavior is known when the classical dynamics is far from 
integrable, namely chaotic (especially, ergodic). The classical ergodic property creates 
a definite semiclassical constraint, expressed by the Schnirelman (family of) theorem(s): 
basically, suitable phase-space measures constructed from the eigenfunctions (e.g., their 
Husimi densities) must tend towards the classical phase-space ergodic measure as ^ 0. 
However, as far as the effective eigenfunction shapes themselves are concerned, this prop- 
erty brings very weak and indirect information, except on the negative side: it is strong 
enough to forbid all familiar semiclassical patterns hke WKB, Gaussian wave packets, 
etc.; but at the same time, it yields no affirmative shape prediction (as yet). Moreover, 
the theorem allows a tiny fraction of eigenstates to evade ergodic behavior, and any of 
these will wholly escape description at this stage. 

We therefore intend to refine the description of quantum eigenfunctions in a classically 
chaotic regime ("chaotic eigenfunctions", for short), and will employ for this purpose 
a holomorphic phase-space representation of quantum mechanics, which can highlight 
certain semiclassical features in great detail. 

We are going to probe a very restricted class of dynamical models which minimize ac- 
cessory complications. At the classical level we take chaotic area-preserving maps acting 
on a 2-d torus (= one degree of freedom), because this is the simplest (compact) phase 
space to admit of such maps with rigorously proven chaotic behavior; out of these we select 
linear automorphisms ("cat" maps) and the "baker's" map. In either case, correspond- 
ing quantum dynamics are readily available. The semiclassical behavior of the resulting 
eigenfunctions can then be examined in this framework, and it already looks quite intri- 
cate. These simple, albeit abstract, models may then lend themselves more readily to 
basic exploratory investigations than more realistic problems plagued with technicalities. 

Our results are mixed, ranging from rigorous to empirical through statistical. We 
definitely do not resolve the internal dynamical structure of eigenfunctions in the chaotic 
regime, but gain several new ideas about their behavior. We now fist these results as we 
outline the contents of the paper. 

The framework just described, and accompanying basic notations, arc summarized 
in Sec. 2, with the main emphasis placed upon the stellar parametrization to be used 
throughout, i.e., the encoding of a wavefunction in the zeros of its Husimi density. In- 
teresting features of the quantum states are not easily extracted from this (nonlinear) 
representation, and some further study of this technique is in fact a parallel topic of this 
paper. 

Sec. 3 reviews the semiclassical ergodic results mentioned above, e.g. like Schnirelman's 
theorem for 'chaotic quantum maps', and lists unsolved questions about the precise way 
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the Husimi densities converge towards the ergodic measure. Various tools (L'^-norms and 
related functionals) can be used to globally measure deviations from strict uniformity, 
and we present extensive numerical results about the Husimi densities of eigenstates for 
the quantum cat and baker maps. In particular we discuss the semiclassical relevance of 
scars in this context. Finally, existing ergodicity-related results on the Husimi zeros (the 
constellations, to be further studied) are recalled here. 

Sec.4 exposes the new rigorous dynamical results of this paper. First, semiclassical 
convergence of the Husimi densities towards the uniform ergodic measure is shown to 
imply the equidistrihution of the corresponding zeros over the phase space in the classical 
limit. We further state as a conjecture that the densities of zeros should do more than 
just equidistribute, i.e., not only do their Fourier coefficients tend to zero, but they should 
decay faster than a certain specified rate. Last but not least, and counter to the belief that 
the phase information of quantum wavefunctions is irretrievable in presence of classical 
chaos, we deduce (from the cquidistribution of zeros) a definite universal behavior for a 
phase-related quantity, the local momentum of the eigenfunction in Bargmann form. This 
result can be cast in a very WKB-like expression for the eigenstate, but to be taken in a 
much weaker, measure-theoretical, sense than for classically integrable systems. 

Sec.5 is devoted to statistical analyses of the distributions of Husimi zeros. The dis- 
cussion continues earlier works based on random ensembles of functions, in the following 
new directions. The results for random polynomials over the Riemann sphere are system- 
atically transported to the torus, for ensembles of theta- functions. At the same time, the 
effects of parity and time-reversal symmetries are separately taken into account. Finally, 
emphasis is put on the resulting statistical properties of the Fourier-transformed densities 
of zeros, with the analytical evaluation of a corresponding form factor (alias structure 
function) . 

In Sec. 6, miscellaneous properties of these Fourier coefficients are observed for eigen- 
states of the quantum cat and baker map, and compared with the dynamical and statisti- 
cal expectations. We find that the previous random form factor imprints each eigenstate, 
through a globally strong suppression of most lower-range Fourier coefficients (up to a 
distance < ^/N). Moreover, we identify a specific signature, within this dual approach, 
of high peaks in Husimi densities (including scars by fixed points of the classical map). 

Sec.7 provides a summary and some overall conclusions. 

2 Quantum models over the torus 

Generally speaking, the most elementary Hamiltonian dynamical systems which can ex- 
hibit chaos are area-preserving maps over a 2-d compact surface. Because such a phase 
space is not of the standard type (i.e., the cotangent bundle of a position space), quantum 
mechanics is not canonically defined; the simplest alternative quantization method needs 
an extra 1-d complex (Kahler) structure chosen on the surface. These compact surfaces 
are topologically classified by their genus g. The most elementary case g — (the sphere) 
is the relevant phase space for spin dynamics [20]. The next case g = 1 (tori) will focus 
all of our attention since it is the simplest to carry proven chaotic maps (the cat and 
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baker maps), while its kinematics (both the classical and the quantal) remain manage- 
able. (The more intricate g' > 2 cases would also deserve investigation, even though the 
physical relevance of such phase spaces is unclear yet; see subsection 2.1.1 for a sketchy 
introduction to quantum mechanics on such surfaces). 

On a 2-d torus, all incquivalent complex structures are labelled by a modular pa- 
rameter, r G {Q'(t) > 0}/SL(2,Z). Here, however, only one r-dependent feature will 
be examined (see Sec. 3.3 and App. A). Otherwise we will only consider the torus 
built from the unit square, with canonically conjugate real variables x = {q,p) G [0, 1)^ 
(mod 1), and fix the complex coordinate z {q — ip)/\/2 e C (corresponding to the 
choice T — i, the most suited to the discrete rotational symmetry of the square; we call 
Tc the resulting complex torus, to differentiate it from the underlying real torus T^). We 
now briefly recall the basic notations and facts pertaining to the quantum mechanics over 
this phase space, fixed once for all. 

Admissible quantum wave functions l'^) are to satisfy two quasiperiodicity relations 
[20, 11]: 

{q+m = e^-'^M?!^) (2) 
(p+l|V^) = e-^i-'^^ (pIV^) 

and they form a Hilbert space TYat ,^ of finite dimension N — 1/h, a condition which 
restricts Planck's constant itself {h = 271%) to inverse integer values; the Floquet angle pair 
[(pi, 1P2) is arbitrary on the dual torus and this freedom makes for a family of quantizations 
which can be labelled by the complex combination (p (fi + 'np2. However, any in- 
dependent effects will be neglected here, being of higher order in h than the leading 
semiclassical features which already challenge our understanding. 



2.1 Holomorphic and phase-space quantum structures over Tc 

We now quickly review (see [20] for details) quantum representations which explicitly 
depend upon the choice of complex structure on T^, here fixed as (the complex torus 
defined above). 

The Hilbert space of eq.(2), Tijv,^ ~ C"^, can also be spanned by doubly periodicized 
Gaussian coherent states, i.e., 

\z)n^^ (-l)^"^'*e2''^('^i"-'^2"*) QiirnQ-nP) (3) 

n,m6Z 

where \z) =^ exp(^a''"/^) |0) is a Weyl coherent state at the point {q^p) of the plane 
(a^ = 2~^/^(g — ip) is the usual creation operator and |0) the ground state of the harmonic 
oscillator); the state 1^;) is thus made to depend holomorphically upon the coordinate 
z = (q + ip) / \/2, instead of being normalized (and {z \ z) = c^^/^). Then for all ip E I^n,^^ 
ip{z) = {z\ip) is an entire function of z, subject to the two quasiperiodicity relations [20]: 
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These functions ip{z) also span an iV- dimensional Hilbert space of entire functions, con- 
stituting the Bargmann representation of TYat,^. Over the torus itself, however, they 
do not define intrinsic functions, only holomorphic sections of the complex line bundle 
specified by the two equations above. It is then worthwhile to seek other representations 
based exclusively on genuine phase-space entities. 

The most popular such construct is the Husimi density, 

H4q,p)=H{x)''^'\{zm'e-^'/'' (5) 

which is strictly doubly-periodic hence defines a (positive) density on the torus; the nor- 
malization condition for the quantum state lip) in 7iN,if translates as 

/ H(x)d^x^l, (6) 

and H{x) can be viewed as a phase-space localization probability for the quantum state 
lip). The Husimi representation thus maps the unit vectors of H.^^^ into probability 
measures over T'^; however, the infinite codimension of this embedding makes it very 
uneconomical, and impractical to invert. 

A contrario [20], we emphasize that ip{z) can be mapped onto a much leaner phase- 
space structure, namely the skeleton of its zeros {z^} (counted with their multiplicities 
and modulo the periods; this is also called the divisor of ip{z)). Any Bargmann function 
has exactly N zeros on the torus Tc, linked by a single linear constraint, 

V2j2^Zk^N^-^j~iip mod[l,i]. (7) 

and conversely, all such A^-uples constitute an explicit and 1-1 parametrization of this 
Bargmann space, thanks to this special case of Hadamard factorization [3, p. 22] for entire 
functions of order < 2, 

iPiz) = {z\ilj)N,^ = const. X e2--2:.=i-'= Uxiz-z^), (8) 

fe=i 

where 

X{z) = 2^^ e'^^' 9i{7rV2z \ i) = (27r)-^/^r(l/4) V(^ | 2'^/^, (9) 

is the Bargmann function of the "singlet" state %, the unique state of Ti.N=i,<p={i+i)/2, 
possessing a single zero at 2; = (we use the notations of [1] for the Jacobi theta functions 
and Weierstrass a function). 

This immediately implies a nontrivial result, in the form of an explicit criterion for a 
Husimi density to be a pure quantum state! Namely, by eqs. (5-8), a Husimi density 
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is constrained first to vanish (quadratically) precisely at the N zeros Zk of ip{z), then to 
satisfy the factorization identity 

N 

H^{x) = n H^{x-Xk) {xk = {qk,Pk) = real form of Zk = 2~^/^(gfe - ipfe)); (10) 

k=l 

> is a constant factor to be adjusted via the normalization eq.(6). Conversely, 
any H{x) having the form (10), for an arbitrary subset of points on the torus, is the 
Husimi density of a wave vector, specified by eq.(8) within the space TCN,ip where now ip 
is constrained by eq.(7). 

A parametrization of pure quantum states is thus provided by this direct 1-1 cor- 
respondence of their Husimi densities with all iV-uples or "constellations" of zeros, or 
equivalently with the (singular, normalized) 2-d Dirac distributions of these zeros, 

N 

pix)"^ N-'Y.^\^-^k). (11) 

k=l 

This new density p{x) can now be studied in place of, or in parallel with, the Husimi 
density H{x) itself. We hope that this nonlinear but "naked" phase-space description 
of pure quantum states, which we call "stellar representation", will shed new light on 
the semiclassical eigenfunction problem. At the same time, however, we must still im- 
prove our understanding and control of this stellar representation which has not yet been 
exhaustively studied. 

The basic relationship between the densities p and H is deduced from the factorization 
property (10), as 

A\ogH{x) = AnN{p{x) - 1). (12) 

This approach thus privileges the logarithm of the Husimi density, as the electric potential 
generated by a point-like unit charge at every zero minus a uniform charge density exactly 
restoring global electric neutrality. 

Now, translating any property of one density concretely in terms of the other becomes 
a central task, but an arduous one. First, the basic relationship (12) is highly nonlinear. 
Then, each zero is a maximally quantal object (contributing one quantum of phase to the 
Bargmann function il^{z)), hence any semiclassical effects (which show up as high-density 
regions) have to be coherent manifestations from a large number of zeros, as observations 
will confirm. 

Due to all the aforementioned difficulties, the stellar representation still requires de- 
velopment; the present work should then be viewed not only in terms of raw results where 
it is perhaps exploratory, but also in terms of methods where it suggests several novel 
approaches and directions for further study. 

2.1.1 Other Riemann surfaces 

As we already remarked, a sphere is the natural phase space for the dynamics of angu- 
lar momentum (or spin). Its quantization involves finite-dimensional representations of 
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SU(2). A stereographic projection of the sphere then allows to represent each state by a 
polynomial (its Bargmann function), leading to both Husimi and stellar representations 
[20]. 

Any higher genus {g > 2) surface Xg can also be given a Kahler structure, as a quotient 
of the Poincare unit disk by some Fuchsian group. Although we are not aware of any 
interesting classical dynamics on such a phase space, we briefly describe the associated 
quantum mechanics to be able to generalize our rigorous results of section 4 to such 
geometries. 

The quantum states on Xg can be defined directly in the 'Bargmann representation': 
they consist of holomorphic sections on some canonically defined line bundles over Xg, 
called automorphic forms; for each allowed value of h, they form a finite-dimensional 
Hilbert space [5], [6, sec. 2.1.2]. Prom these holomorphic sections, one can also define 
Husimi densities, which verify factorization properties similar to eq.(lO), so that a quan- 
tum state will be defined uniquely (up to a global prefactor) by its constellation of zeros 
on Xg [6, sec. 4.6]. 

2.2 Semiclassical eigenfunction problem in reduced setting 

The issue stated in the introduction can now be considered in the special setting of 
quantum maps on the torus. Let us mention that similar considerations work on the 
2-sphere, which is the phase space of a classical spin [20]. 

We first rapidly comment on the problem of quantizing a given area-preserving map 
on the torus, i.e. how to associate to this classical transformation a sequence of N x N 
unitary matrices Un,,p acting respectively on Hilbert spaces Ti-N.ip, with correct semi- 
classical properties [7, 8]; the whole sequence of matrices {Un,:p}n>i is called a quantum 
map. 

Such a quantization is in general not uniquely defined, assuming there exists one. A 
systematic way of quantizing symplectic maps was recently devised by Zelditch, using 
a Toeplitz operator formalism [9]; however, this formalism does not in general lead to 
closed-form expressions for the Un,,p, and is currently restricted to smooth maps. On the 
other hand, if the classical map is piecewise affine, a 'down-to-earth' construction provides 
explicitly a sequence of unitary matrices. We will restrict our numerical investigations to 
maps of the latter type, namely the baker's and the cat maps, which are already quantized 
and studied in the literature [38, 39]. 

Whereas the Ploquet angles (/? will be fixed and subsequently implied, dependence upon 
the integer N will be a crucial feature for us, with the semiclassical regime corresponding 
to = 1/h +00. 

(Prom the standpoint of physical systems, which are described by quantum Hamilto- 
nians as opposed to maps, the above construction crudely models a quantum-dynamical 
reduction of a 2-d time-independent Schrodinger equation to a 1-d quantum map, analo- 
gous to a Poincare surface-of-section map for a classical fiow [35, 36]; our problem is thus 
'of reduced type'.) 

The issue is then to compute (or describe) all (or some) of the eigenfunctions of the ma- 
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trix family {Un} asymptotically as iV — > +00, in connection with the classical dynamics. 
We will moreover focus on phase-space descriptions for these wave functions, privileging 
the stellar representation for its nonredundancy; we coin the shortened notations "eigen- 
density", resp. "eigenconstellation" , for the Husimi density of an eigenf unction, resp. its 
constellation of zeros. 

The two quantum maps we select both correspond to fully chaotic area-preserving 
transformations of the torus in the classical setting. One is the baker's map, quantized 
for even N [39, 40], and under antiperiodic boundary conditions {(p = (l-|-i)/2 in eq.(2)) so 
as to preserve the classical parity symmetry P: {q 1— > —q, p 1— > —p}. The other consists of 



(with integer entries and unit determinant) such that tr(5') > 2; these will be quantized 
under periodic boundary conditions (i.e., (p — 0) [38]. Both of the resulting quantum 
models commute with spatial parity, hence their eigendata separate into even and odd 

subsets. 

As concerns 'cat' maps, moreover, we will look at their quantizations only for certain 
prime values of N that are 'splitting' [42] and such that spectral degeneracy (a potential 
source of extra complications) affects just one eigenvalue (of multiplicity 2); for these 
values, semiclassical ergodic behavior is fully proven [12]. Wc will however have to switch 
between two distinct classical cat dynamics, of respective matrices 



The reason is that S gives a time-reversal-invariant dynamics (a consequence of a = 
d in that matrix), and this quantum-mcchanically implies real eigenf unctions. While 
this symmetry (T) does not interfere with most computations, it greatly complicates the 
analysis of random ensembles to be used as comparison models (Sec. 5); at that point 
it becomes simpler to work with the other map like S', which has no anti-canonical 
symmetry and is as easy to quantize, using the formulas in [38] . To prove that S' has no 
extra symmetry (besides P, always present), we had to study conjugacy classes in SL(2, Z), 
and found no such map with a trace < 36. More generally, in order to avoid anti-canonical 
symmetries, in practice one has to accept much larger Lyapunov exponents, which in turn 
make numerical studies of semiclassical behavior (e.g., the search for scars) much more 
delicate. Consequently, the map S will be more suitable for asymptotic verifications, and 
S' for statistical ones. 

3 Measures of semiclassical ergodicity 
3.1 The Schnirelman property 

Semiclassical ergodicity reads as follows for quantum maps on the torus. 

Let {^f )}, 

=i....^N be orthonormal eigenfunctions of a quantum map 11^- defined over 
the space Hn- Let C denote the (normalized) Liouville density over the phase space for 



the cat maps or linear hyperbolic automorphisms, given by matrices 





(13) 
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the corresponding classical map (namely, the invariant area density) . If this classical map 
is ergodic with respect to then almost any subsequence of Husimi "eigendensities" 
{H (N) } converges to C as N ^ +00, in the weak-* topology for measures. 

(Weak-* convergence of a sequence of measures {/in} on a compact space X is defined 
by the convergence of every scalar sequence fx f /-I'M where / is any Gxed continuous 
function [45, p. 113].) 

We will say that a sequence of measures {/xat} ^* C "has the Schnirclman property" — 
by reference to the name "Schnirelman's theorem" generically used for the many dynam- 
ical versions of the whole statement above (the first one was worded for Riemannian 
Laplacians [10], while we directly stated the quantum-map version). 

In the theorem, "almost any" means that the stated property may be violated by an 
exceptional subset of eigenfunctions of asymptotic density zero, i.e., if this subset reads 
as {^f ^},eo(iv) then #{l](iV)}/iV ^ 0. 

For the torus maps we work with, the invariant measure C is just the usual flat measure 
of uniform density 1 on . Then the Schnirelman property for a family of Husimi densities 
is their uniform spread (as measures) over the torus as ^ +00. 

In full mathematical rigor, the above theorem was proved in [9] for arbitrary smooth 
ergodic classical maps (on the torus), quantized through the Toeplitz formalism; this 
includes in particular quantum cat maps [11]. Moreover, direct calculations establish that 
the particular sequences of quantum cat eigenstates we will numerically study have the 
Schnirelman property [12]. By contrast, in the case of the baker's map the theorem is 
only partially proved [13], but is believed to hold nevertheless. 

We now display (fig.l) a sample of Husimi eigendensities of quantum cat and baker's 
maps. Already at a glance, the situation shows a wealth of complex behaviors which 
cannot be exhausted by just invoking the Schnirelman property. The densities do show an 
overall tendency to uniformization but only in a very coarse sense and at quite slow rates. 
Subdominant structures are quite significant, in fact their intricacy and profuscncss wildly 
grow with N. Certain patterns seem to occur recurrently but with no definite regularities; 
in particular, periodic orbits are often scarred but not in an obviously systematic way. 

Our problem of concern will be, in general terms, to find various ways to tackle on a 
finer scale such eigenstates obeying the Schnirelman theorem (which constrains them on 
the coarsest scale only). This includes, e.g., the popular topic of scarred eigenfunctions: 
how they should be derived, and how they fit with the Schnirelman picture. Unresolved 
issues in this area are: 

1) sharpen the Schnirelman property: 

la) assess convergence in stronger topologies; 

lb) seek equivalent conditions upon the eigenfunctions themselves; 

2) describe the finite- iV situation: 
2a) estimate the convergence rates; 

2b) analyze and classify the quantum corrections H jn) — C and in particular any excep- 
tional, non-ergodic ones. 
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3.2 Norms and other invariant measures of density fluctuations 

One natural way to estimate the fluctuations of a density H{x) away from flatness is by 
L''-norms and related functionals. Keeping in mind that H[x) is a nonnegative density 
on the torus (and that the latter has unit area), one defines, following [16] 



\u\\ '^^^ 

\-ti\\r — 



d'xHixY 

T2 



l/r 

(r > 0). (14) 



For r > 1 these quantities define norms, called U , out of which the L^-norm is taken 
by eq.(6) to fix the normalization, as \ \H\\i = 1; then by convexity (Jensen's inequality, 
see for instance [44, p.70]), \\H\\s > \\H\\j. > 1 for s > r > 1, with equality precisely 
characterizing the fiat case H{x) = 1. Finally, those (and derived) quantities are geometric 
or shape invariants. 

Intuitively, fluctuations of H{x) push ||-ff||r upwards, more so when r is larger. Ulti- 
mately, 1 1 -f^ I loo (= -f^max, the sup-norm) records the highest enhancement only. A more 
global assessment of the non-uniformity requires an intermediate norm, of which the sim- 
plest is the L^-norm; since H is quadratic in the wave vector, is an integral quartic 
in ip, like the inverse participation ratio commonly used as a measure of quantum local- 
ization (but often in conflguration space). 

The weakest norm, hence "least remote" from the weak-* topology, would be the L^- 
norm, had it not already been used up for normalization. The closest approach then 
considers \\H\\r as r ^ 1+, or rather the functional (non-negative, by convexity) 



"^llffll' 




^ \\r 

dr 


r=l 



d xH{x) logH{x) 



which has another important meaning, as the information content carried by the probabil- 
ity density H (Boltzmann's H-function, or negative entropy). This quantity is related to 
the 'classical-like' entropy of a quantum pure state (more generally, of a quantum density 
operator) introduced by Wehrl [17], then studied (on the plane and the 2-sphere) by Lieb 
and Lee [18]. This function can be used to deflne the 'dynamical quantum entropy' of 
a unitary map [19], which is supposed to generalize the Kolmogorov-Sinai entropy of a 
classical system. 

Last but not least, since the Husimi functions of pure quantum states (the only ones 
of interest here) have a multiplicative structure, the geometric mean of the density H is 
a natural functional to use [44, p. 70], 

GM[//] = exp/ d'^x\ogH(x) Mm \\H\\r) (15) 

Jt2 r^0+ 

since it is itself multiplicative {GM[HiH2] = GM[Hi]GM[H2]) ■ Thus, for instance, the 
evaluation of GM on both sides of the Husimi factorization formula eq.(lO) fixes the 
normalization factor C^, resulting in a more precise expression for this factorization as 

AT 

H^x) = GM[H^] n H^{x - Xk); H^x) = H^{x)/GM[H^] (16) 

k=l 
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(with GM[H^] = (27r) -3/21(1/4)2; see App.A.l and table 1 below). More generally, a 
Husimi density will be denoted H when we normalize it by its geometric mean (instead 
of L^-norm). 

GM[i/], like ||-ff||r for r < 1, is not a norm (all convexity-related inequalities get 
inverted with respect to r > 1); it satisfies < GM[i7] < 1, with equality in the flat case 
whereas fluctuations of H drive it downwards. 

In summary, the L^- and sup- norms, the H-function, and the geometric mean, all give 
useful size estimates for the fluctuations displayed by flnite-A^ Husimi densities. Although 
they yield largely similar density rankings, the various functionals are inequivalent; at the 
same time we see no compelling argument to favor any single one of them. 

3.3 Comparison states 

As we intend to probe the finer structure of Husimi eigendensities, beyond the Schnirelman 
property H ^^-4 and since those are pure states, we simply cannot dismiss their mul- 
tiplicative structure at the finest quantum scale: any such Husimi density is the product 
of N translates of a single elementary building block, and is thereby determined by its 
constellation of zeros (up to normalization: cf. eq.(16)). Thus, quantum mechan- 
ics imposes a granular microstructure upon these densities at finite N, forcing them to 
vanish at precisely points. If pure-state Husimi densities uniformize as N oo, it 
can then happen in a coarse-grained sense only; in no way can they look smoother and 
smoother, "classical- like" , in the semiclassical limit, as their inhomogeneities will only ex- 
perience a shrinkage of their geometric (x-) scale and not of their intensity. Any fine-scale 
semiclassical theory must then heavily acknowledge such rigid structural constraints. 

A first step in this direction is to gather data about various patterns of zeros that 
are at the same time characteristic enough and analytically tractable. Such reference 
constellations will then provide a comparison scale: we proceed to define a few of them 
now (up to unspecified translations). 

3.3.1 Lattice states 

Intuitively speaking, the most equidistributed Husimi densities subject to the constraints 
(16) must be found among lattice states, whose zeros form a 2-d lattice A with a funda- 
mental cell of area 1/A^. Such displays of zeros have maximal, "solid-state" rigidity. Some 
of them occur in exceptional eigenstates of cat maps. 

Up to scale invariance, the shape of a lattice A can be described by a complex number 
T with $J(t) > (the modulus) such that the torus C/A is conformally equivalent to the 
complex torus T^- of generators {1, t}. A Husimi density specified by the set of zeros A can 
be directly written in terms of the (properly rescaled and displaced) elementary (A^ = 1) 
Bargmann function xi^ I t) quasiperiodic w.r.t. and its Husimi density i7^(^)(g,p), 
which are expressed by the formulae [54] 

Xiz I r) = (2$J(r))^/^ e^ ^i(7rV2^ | r), (17) 
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ip) I r) 



the normalization being 




area(ri-)). 



(18) 



Although the present article restricts to quantum mechanics on the square torus 
Tc = Ti (for which we recover our standard xi^) ^-iid H^{x)), this extension of scope 
to functions quasiperiodic w.r.t. any torus T^- is useful for three reasons. Firstly, Husimi 
densities of the form H^(^T--j{Nq, Np) for rational values r = {m + in)/N have the square 
torus periodicity too, and lattice eigenstates of quantum cat maps have precisely this form 
(section 4.4). Secondly, in section 4.2.4 we will define Husimi densities whose constella- 
tions are deformed lattices, and their asymptotic expressions will make use of the function 
Hx{r){(l,p) for variable r. A third reason (already stated but excluded from the scope of 
the present work) is that varying r (over the modular domain) also amounts to exploring 
the range of complex structures over a fixed real torus, or all possible Bargmann-Husimi 
representations for a prescribed problem. 

For lattice states, the various norms and functionals are functions of r alone (by scale 
invariance). With the phase space itself being of square shape, the invariants for r = i 
arise the most readily (they describe the elementary Husimi factor of eq.(lO): see 
figs. 2a, left and 4, left); however, we expect the optimal equidistribution to be attained 
by the closest packing of the zeros overall, namely the equilateral triangular lattice (i.e., 
r = e''^/^). While this geometry cannot be realized on a square torus, it can be approached 
in the limit N — > +00 (scale-equivalent to the infinite-plane limit; cf. fig. 2a, middle), so 
we conjecture that its invariants will be the good ultimate bounds. 

(In fact, we seem to encounter natural extensions of the close-packing class of problems: 
to find the constellation(s) of N points on the torus that optimize any one of the invariants; 
the difficulty may increase by changing the surface, as the review [49] on the spherical 
case suggests.) 

We can already prove (see Appendix A) that two of the above invariants (namely, the 
geometric mean and the L^-norm) are optimized by the equilateral lattice within all lattice 
states. Curiously, such invariants are related to spectral determinants of tori, studied for 
instance in [50]. 

3.3.2 Localized states 

At the other end, we should consider patterns with highly concentrated zeros. Examples 
forming an important class are the eigendensities of classically integrable systems, for 
which both the Husimi densities and the densities of zeros respectively concentrate on 
lines: one prototype is a plane wave (or delta wave, up to 7r/2 rotations), which is in fact 
just a highly anisotropic lattice state (having r = i/N: fig. 2b, left). Still more localized 
is a coherent wave vector defined by eq.(3), with its Husimi density concentrated around 
a point in phase space; its zeros lie on two perpendicular axes (those which meet at the 
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antipodal point: fig.2b, middle). The most extreme concentration is however realized 
by the completely degenerate pattern in which all N zeros coincide (fig. 2b, right); even 
though this state is mainly a curiosity and has a Husimi density barely different from the 
previous one, we include it for the sake of completeness. 

3.3.3 Statistical states 

In-between the previous two classes, we may consider statistical distributions of zeros for a 
random polynomial ensemble, inasmuch as these reproduce quite well (albeit empirically) 
the local fluctuations seen for zeros in eigenfunctions of classically chaotic maps. As in 
the random-matrix modeling of spectra, the dynamical symmetry of the problem must be 
taken into account. The random coefficients have to be taken fully complex for a model 
without any antiunitary symmetry (a sample on fig. 2a, right), and real otherwise (see Sec. 
5 for details). 

3.3.4 Analytical results 

The norms and invariants of the previous comparison states can be computed: either 
exactly, or asymptotically for large N, or just numerically for a few. The results are listed 
in table 1 (upon setting all L^-norms to unity, as usual). 

The lattice-state invariants arc computed long-hand using classical theta-function 
identities (see also App. A). A tricky point for the sup- norm is that the maxima of 
the Husimi density H^(^t-) depart from the dual lattice when r ^ i (for r = e''^/^, they lie 
at the centers of the equilateral triangles of zeros). 

The concentrated state examples (last three fines) are also expressed using theta- 
functions, and furthermore all required integrations are performed asymptotically (by 
stationary phase). Remark: the study of the Bargmann representation over the whole 
plane (i.e., for the Hilbcrt space L^(k)) shows [18, Theorem 3] that for each r > 1, the 
L^-norm of is well-defined for all ip G L^(m) and attains its maximal value when ip is 
a Weyl coherent state (and fikewise for the H- function). It is reasonable to expect these 
results to extend to the torus, at least semi-classically. Then, with our normalization, 
those invariants would be bounded as follows over Hn'- 

Vl<r<+oo, sup \\H^\\/J^ Mr{N) -JL^ {N ^ oo) (19) 

||V||i=i (iVrjV'- 

sup R[H^] = Mi{N) ~ logA^ - 1 (A^ ^ oo) 

1^11=1 

and these maxima Mj.{N) should be reached at (or near) a torus coherent state; indeed, 
their asymptotic estimates agree with those of the corresponding coherent-state invariants 
in table 1. 

The random-model averages are computed from the distribution law for the Husimi 
density values, P{H) ~ in the N — > -|-oo limit (for either the real or the complex 
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state 



sup-norm L^-norm 



H-function Geometric Mean 



equilateral 

lattice 1.4610 ^ 1.0768 0.1016 ^ 0.8435 

SQuare '^(V^)' r(i/4) r(i/4)^ 

square V2 7r3/2 ^2^3/4 (27r)3/2 

lattice fti 1.6693 ^ 1.0864 0.1110 r; 0.8346 



random <21ogAr (1 - 7) ± V^^^^^^4^ e-(l ± ^M) 

plane wave ~ V2iV ~ AT^^ ~ i[log(2iV) - 1] ~ ^2iVe-'^^/6 

coherent ^ N ~ ~ log - 1. ~ Ne'^^/^ 

degenerate ^ N ~ ^iV/2 log - 1 ~ TV 2"^ 



Table 1: Values of the various fluctuation measures for several comparison states, from 
the most ergodic (first two rows) to the most localized (last three rows). In the third row, 
we give typical values for a random state (7 is Eulcr's constant), as well as the standard 
deviations (numerical values of the averages appear on fig. 3a). 

Gaussian ensemble). We can then compute any integral in average as 

([ f{H{x))d'x) ^r^f{H)e-^dH (20) 
\jt2 In Jo 

and this yields the constant values listed as typical for the ensemble. The standard 
deviations from the average values are estimated using integrals of the type: 




for real positive {s,t) (see Sec. 5 for details). The sup- norm expectation value diverges 
in this regime, and needs a more delicate finite- analysis (App. C). 

The "most naive" ergodic assumption (i? ~ 1 pointwise) would set all invariants to 
unity (except the H-function, to zero). By contrast, the tabulated values for the lattice, 
resp. random, models mark the incompressible, resp. the typical, levels of quantum in- 
homogeneity. Far from vanishing in the semi classical limit, these levels are independent 
of A^: this rules out convergence towards classical ergodicity in any of the corresponding 
stronger topologies. Furthermore, the random invariants are substantially above the op- 
timal lattice values, confirming the visual impression from the two rightmost plots on fig. 
2a that random densities are quite non-uniform. 



3.4 Invariant data for eigenfunctions 

Concerning now the invariants of eigendcnsitics corresponding to classically chaotic dy- 
namics, our results consist of numerical data displayed graphically. We have computed 
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the full set of invariants for all eigenstates, at selected values of N, of two quantum maps 
with highly chaotic classical limits, i.e., the baker's map and the cat map S specified in 
Sect. 2.2. We plot and compare the distributions of the invariants for the resulting eigen- 
state families. Fig. 3a plots the values of four functionals for the even-parity eigenstates 
of the quantization of the cat map S, with the geometric mean plotted downwards for full 
visual consistency (the odd-parity plots look much the same, and are not shown); Fig. 
3b plots the values of two functionals for both parity eigenstates of the quantum baker's 
map. 

3.4.1 General observations 

The figures first show the values of the invariants to be substantially depressed in com- 
parison to the concentrated state examples, but still broadly scattered, with relatively 
few isolated values or clusters: this confirms the visual impression that individual eigen- 
functions within a single matrix (a quantum map at given N) retain widely differing 
global fluctuation patterns. For each eigendensity the value of a particular invariant is 
determined by its N zeros, but large values of this number do not suffice to enforce a 
statistical averaging towards a universal behavior in this respect. Besides, that ordering 
of eigendensities according to their global contrast is generally barely sensitive to the par- 
ticular choice of functional. At a more detailed level: the cat fluctuation values are much 
more squeezed than the odd baker values, themselves more squeezed than the even baker 
ones. This difference could be related to the singular nature of the baker dynamics, whose 
discontinuity will be more felt by even-parity states (especially those scarred above the 
irregular fixed point (0,0) which straddles the discontinuity, whereas the period-2 point is 
regular) . 

We next see that, upon further averaging over the entire eigenfunction set (for a given 
N) , each invariant then comes fairly close to the typical random value (third line of table 
1). This convergence is striking for the cat map, whose eigcnfunctions thus appear to 
fulfill generic expectations (of compliance to a random model) somewhat better than the 
baker's map in spite of its arithmetic idiosyncrasies. (This is the exact opposite of the 
situation for the eigenvalues: cat map spectra are highly non-generic, whereas the baker 
map spectrum is fairly GOE.) The reason may again be that the baker eigendensities 
tend to be unusually bumpy above the map discontinuities. By contrast, the cat map 
possesses a few specially regular lattice and crystal states (always yielding values of the 
functionals among the lowest-lying), but these atypical states do not seem to drive the 
averaged behavior away from genericity (at least for our selected values of A^) . 

We stress again that those density fiuctuations (of the same order as in the random en- 
semble) already represent an appreciable and invariant (A^-independcnt) degree of spatial 
inhomogeneity. We can raise to a conjecture the empirical observation that in average, 
the invariants appear to be reliably constrained to the universal random model values. 
This entails that the average convergence of eigendensities towards uniformity is truly 
poor. Eigendensities displaying a reasonably high level of homogeneity or rigidity must 
then have invariants well below that average (for example, the lattice eigenfunctions in 
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cat maps), and for the sake of global balance some other states will then have to display 
higher inhomogeneity or localization than if they were random (such as numerous states 
scarred by periodic orbits). 

Concerning A^-rclatcd behaviors: most fluctuation sizes grow fairly slowly with A^, 
especially for the cat map; A^-dcpcndcnces are erratic as a rule; their irregularity is reduced 
but not wholly suppressed when a restriction is made to special values of (primes, 
powers of 2, etc.) and/or to specially selected states (most scarred, or most ergodic, etc.). 

All this confirms that the individual behaviors of eigendensities may not be easy to 
characterize or classify, whereas the quantum fluctuations seem to follow the random 
model not only microlocally, but also globally once averaged upon all eigenstates. 

3.4.2 About point scars in eigenfunctions 

The plots allow us to start a quantitative analysis of scars in eigendensities. The sup- 
norm, and to a lesser extent the L^-norm, are the best invariants to select the states 
with the strongest local density enhancements. One currently debated issue concerns the 
semiclassical weight of scars: in particular, do scarred states still follow the Schnirelman 
property or not? This question is hard to settle partly because this convergence property 
is specifically weak-* (and not metric, for instance), hence it is difficult to establish a 
single and intrinsic general scar strength scale. 

We can however get a clearer idea about the simplest, point-like scars. Scars are 
features localized near unstable periodic orbits, and in reduced dynamics these consist of 
discrete points. Quite a few Husimi eigendensities of maps do present sharp enhancements 
that are very localized near periodic points, and we can focus upon the most strongly 
scarred states of this type. We always find the widths of such point scars to scale, very 
rehably, like N~^/'^ in all phase space directions (matching the distances of the nearest 
zeros when these distribute uniformly, as the figures confirm). The value of the Husimi 
density at the scarring point (the sup-norm if the highest scar is under scrutiny, as will 
be the case here) is then a good strength indicator: the area of such a scar being of order 
l/N, a safe criterion for it to produce an asymptotically vanishing contribution (as a 
measure) in a sequence of densities {Hj^} is 

a = limsup(log I li^Arl loo/ log A'") < 1 

W— »oo 

hence | |i/Ar| |oo = o{N): the height of the scar should grow strictly less than the sup-norm 
of a coherent state (Table 1). Observing each log- log plot of sup- norms, we can then 
compare the growth trend (with N) of the uppermost points in the cloud of sup-norm 
values against the (unit) slope of the straight line of coherent state sup-norms, to get a 
(necessarily rough) idea about the asymptotic persistence of such scars. 

For the cat map, those uppermost points are irregularly distributed with respect to N, 
nevertheless heuristic upper bounds for the slope a are 0.5 for odd states (not shown) and 
0.34 for even states; this is consistent with the proven fact that all cat map eigenstates 
used here satisfy the Schnirelman property. For the baker's map, the most strongly scarred 
states form more regular families and estimates look more robust; they yield slopes below 
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0.85 for odd states and 0.8 for even ones (but the even-state scars are stronger in the 
range of N under display). The slopes for the baker's map are thus much closer to the 
critical value a = 1, while wc cannot control the error range of our crude estimation 
process. Hence all we say is that the plots slightly favor the hypothesis that even the 
most strongly point-scarred baker eigenstates still satisfy the Schnirelman property. 

We must however add that there also exist eigenstates which are strongly enhanced 
along whole stretches of stable and unstable manifolds (or neighboring arcs of hyperbolae) ; 
the baker's map shows many such instances, especially about the (singular) fixed point 
(0,0) (fig. lb, right). Because of the extended character of these scars, their sup-norms 
need not lie anywhere as high as those of point scars for them to have stronger classical 
imprints; at the same time, the functionals used above may not suffice to quantify this 
strength precisely enough. 

3.5 Ergodicity and the zeros 

A deeper way to incorporate the multiplicative structure (10) of eigendensities into their 
analysis is however to handle the eigenfunction problem entirely within the stellar rep- 
resentation, or to describe directly the phase-space densities of zeros p{x) of individual 
eigenfunctions. In particular, the semiclassical eigenfunction problem is then to unravel 
the large- behaviors of these "eigenconstellations" . 

Ideally, we would like to be able to to compute or control p{x) (the zeros) ab initio, 
directly from the underlying quantum-dynamical equations, and we have begun investi- 
gating this approach as well. However, the simplest Schrodinger operator already yields 
for the zeros a dynamical system of nonlinear equations involving strong M-body interac- 
tions [22], where M ranges from 1 to min(degree of potential, A^) ; apart from a few trivial 
cases (always yielding completely rigid motions for the constellations), it is currently un- 
clear how such a system can be handled. In our reduced setting, moreover, the dynamics 
is that of a quantum map, i.e., it is no longer local, meaning that zeros must fare even 
worse; and indeed, we cannot in the least simplify the equations of discrete-time motion 
in terms of the zeros (except for certain 'lattice states' under cat maps). 

That leaves us with indirect approaches where we try to relate the stellar representation 
to more familiar ones, especially Husimi's. Similar questions can be asked about p{x) as 
were listed for Husimi densities H{x) in Sec.3.1; now, moreover, any answers will directly 
refer to the eigenfunctions via the stellar parametrization, implicitly settling question 
lb). However, each property now requires translation from one density to the other: in 
particular, the Schnirelman property, the scar phenomenon etc., should be redescribed in 
terms of the zeros, with the hope of a better understanding later. Concrete problems of 
this sort are, typically: 

Ic) describe the densities of zeros and their higher-order correlations, with their de- 
pendence on classical dynamics; 

2c) classify the fiuctuations of their distributions as dynamical (e.g., scars?) or statis- 
tical (and: universal?). 
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3.5.1 Existing observations 

Some features of eigenconstellations have already received dynamical interpretations [20]- 
[24]. 

For integrable systems, eigenfunctions follow a WKB-type Ansatz like (1) in the 
Bargmann representation too, from which it follows that eigenconstellations have to co- 
alesce as ^ oo onto fixed curves, namely certain anti-Stokes lines of the complex 
classical action in the z variable, along which the zeros moreover distribute regularly with 
spacings of order 1/iV. 

In chaotic systems, by contrast, numerical computations of Husimi eigendensities in- 
dicate that eigenconstellations roughly equidistribute in (almost) all of phase space, with 
spacings of order l/\fN] moreover, higher-order correlations between the zeros are found 
to very accurately follow a universal model, namely the statistics of the zeros of a Gaus- 
sian random ensemble of polynomials (which acts nicely as an analog for eigenfunctions 
of random-matrix ensembles for eigenvalues). All those findings about chaotic systems 
have remained mostly empirical. 

More systematic observations of eigenconstellations at increasing A^ (fig-1) show a 
proliferation of vaguely lookalikc but not uniform patterns in want of ordering. It is 
specially challenging to try to correlate some spectacular fluctuations of eigendensities, 
like scars above unstable periodic points, with specific patterns of zeros. Around every 
high peak of the Husimi density there is necessarily a region devoid of zeros, and sometimes 
the nearest of these vaguely affect a hyperbolic distribution, but the connection is tenuous: 
the decrease of the density of zeros under a scar is quite moderate and, conversely, regions 
with an equally low density of zeros do not build high Husimi densities above them as a 
rule. 

3.5.2 Some directions of analysis 

A more systematic description of the eigenconstellations on the one hand, of the high- 
density fluctuations (e.g., scars) within the Husimi densities on the other hand, and 
of explicit relationships between the two, is still lacking. We just know as a matter of 
principle that the zeros have to embody a full description of the states, which then includes 
their semiclassical behavior. The issue is then to locate and decipher the dynamical 
information {order parameters) buried in the zeros. 

At best, each zero in an eigenconstellation might follow some predictable rules, at 
least scmiclassically; but even if wc restrict to exceptional subsequences of states or zeros 
this may be too demanding, as the dynamics of the zeros is exceedingly complicated. 
At worst, the information could be so scrambled up that only quasi-random features 
could be asserted. We will seek to identify intermediate types of information within the 
zeros, weaker and thereby hopefully more accessible than the individual locations, but still 
dynamical (i.e., stronger than statistical). These could lie within collective coordinates of 
some sort (e.g., we recall that semiclassical effects can only be collectively generated). 

One possibility inspired by the integrable situation (regularly spaced zeros on a curve 
produce WKB-like wave functions) would be to build up wave functions through patching 



18 



together local 2-d patterns of zeros of a few definite types. This approach is valuable if 
forms of short-range order can be identified in eigenconstellations even when the classical 
dynamics is chaotic. Such regularities are actually not infrequent, but are not systematic 
either: chaotic eigenconstellations mostly look like gases or liquids of zeros, as opposed 
to solids (except for the few lattice eigenstates in cat maps). Zeros thus appear to sub- 
stantially use their 2-d freedom of relative motion, making it difficult to greatly reduce 
the number of their parameters. 

If individual zeros prove difficult to isolate, a valuable complementary approach is to 
analyze their density p{x) in the dual Fourier space which, for our torus phase space, is 
the lattice of integer vectors k = {kg, kp) e Z^: the Fourier coefficients of the density are 

1 ^ 

Pk = i^Y.^~'^''-"' ^ez2, (21) 

under the normalization pq = 1. Some valuable features of this approach are, for instance: 

1) the basic relationship (12), from the logarithm of the Husimi function to the density 
of zeros, becomes diagonal in Fourier space: denoting 

h{x) = logH{x) = J2 hk^i^^'^-'' (with /lo = 0), (22) 

we find that eq.(12) is mapped to the simple relation 

Pk-k = -T:{kl + kl)hk^-T:\k\^hk {yk&I?), (23) 
where 5^ 1 if A; = (0,0), else (Kronecker symbol); 

2) if the zeros form a lattice, then the Fourier transform is simply the Dirac delta 
distribution on the dual lattice (this expresses the Poisson summation formula); 

3) if the zeros form a more disordered distribution of points, this Fourier transformation 
can be well controlled statistically in terms of structure functions, alias form factors (see 
Sec. 5). 

However, this Fourier analysis is deferred to Sec. 6 of this paper, because some criti- 
cally relevant information still has to be gathered before: some rigorous dynamical results 
about the zeros are presented next, whereas a statistical approach is developed in the 
Section thereafter. 

4 Dynamical results 

In this section we study the analytical properties of sequences of states with the Schnirelman 
property, i.e. such that their Husimi densities converge weak-* towards the Lebesgue 
measure C on the torus. We prove that the constellations also equidistribute in the weak-* 
sense, and we deduce from this some universal properties of the phases of their Bargmann 
functions. We also give examples showing the limitations of this rigorous approach, as well 
as stronger results for a class of eigenstates of quantum cat maps (i.e. lattice eigenstates). 
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4.1 Equidistribution of the eigenconstellations 

Using the fact that the Husimi density is the square of a holomorphic function (up to a 
trivial Gaussian factor), it is possible to extract valuable information about the densities 
of zeros Pn{x) from Schnirelman's property. 

The link between both densities is provided by the logarithm of the Husimi density. 
Namely, the factorization property (16) implies the very simple relation 

N 

\ogH{x) =J2h{x- Xj) = N{h * p){x), (24) 

where the function h{x) is the logarithm of the building block H^, shifted such that its 
average over vanishes: h{x) = \ogH^{x). This function corresponds (up to the sign) 
to the electric potential generated by a delta-like unit charge at the origin, balanced by a 
uniform negative charge to ensure global electric neutrality. Precisely, we have 

Ah{x) = 47r{6{x) - 1), (25) 

from which we easily recover eq.(12). 

Using the above representation of the Husimi functions, we will now prove that 
Schnirelman's property for a sequence of Husimi densities implies the weak-* convergence 
of the corresponding densities of zeros to the Lebesgue measure £ as A?" — > oo. 

The densities of zeros are positive and normalized, so they belong to a compact set 
of measures in the weak-* topology (this is the Banach-Alaoglu theorem [45, p. 105]). 
Therefore, we can extract a subsequence converging to a positive normalized measure Poo 
(the elements of this subsequence will henceforth be noted p^v)- 

Our first task is to compare as precisely as possible the finite- A?" and limiting electric 
potentials Hn = h * pN and hoo = h * poo- The weak-* property pjv — > Poo obviously 
implies h]\f ^* hoo, but we need information in stronger topologies. 

h{x) is at the same time upper semi-continuous (u.s.c) [44, p. 37] and in L*(T^) for 
all 1 < s < cxD. This, combined with the fact that the measures p# are normalized and 
positive (# stands for N or oo), implies that the potentials are also u.s.c. and in 

L^(t'^), with < \\h\\s [46]; furthermore, they are bounded above by M =^ suph{x) 

(all these properties can be shown by considering a decreasing sequence of continuous 
functions converging pointwise to h{x), as in [46, theorem 3.6]). 

These properties, combined with the weak-* convergence of the pjv, imply [46, p. 2 09]: 

Vx, hoo(x) =^ limsup /iAr(a;) exists, 

and hoo{x) < hoo{x). (26) 

Equivalently, Vx, h^ix) is smaller than hoo{x) + e for large enough. This inequality can 
be proven to be uniform w.r.t. x by using subharmonicity properties of the potentials. 
The functions themselves are not subharmonic, because their Laplacians involve a 
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uniform negative charge (independent of N). Removing this constant charge amounts to 
systematically adding the function ttx^ to the potentials. The resulting (non-periodic) 
potentials g^{x) — h#{x) +7tx'^ are uniformly bounded above over the square [0, 1]^ (like 
h^), but they are also subharmonic. Therefore, we can apply to the sequence {qn} several 
lemmas pertaining to subharmonic functions (see for instance [47, §1.8] or [48, chapter 
3]), and then transfer them back to the /i-potentials by subtracting the ttx"^ term. 

To obtain the uniformity in eq. (26), we need to consider a decreasing sequence of 
continuous functions {/ioo,m}meN on converging pointwise to h^{x) (the existence of 
such a sequence is due to the upper semi-continuity of h^). Then, we have Hartogs' 
lemma [47, theorem 1.31]: 

Lemma 1 Vm e N, Ve > 0, 3N{m,e) s.t. 

N > N(m, e) =^ Vx e T^, hN(x) < hoo,m{x) + e. 

Before using this uniformity result, we state additional properties of the function 
hoo{x)i also due to the subharmonicity and the uniform upper-boundedness of the g- 
potentials. 

Lemma 2 The function x i— > g^{x) = hoo{x) + ttx"^ is "almost subharmonic" . 

This means that it is equal almost everywhere to a (unique) subharmonic function g^^. 

This function can be defined as an upper regularization ofg^: glo{x) — lim^-^o I{\y-x\<r} ^oo(y)^Z/) 

which yields the property ^oo(^) — ^oo(^) 

gl^ can also be defined as the lowest subharmonic major ant of g^. Therefore, from 
the inequality (26) we deduce 

hoo<Ko< hoc. (27) 

This, combined with the weak-* convergence ^ h^, entails the identity h*^ — h^o- 
Therefore, (26) is actually an equality for almost all x. 

In a second step, we use lemma 1 to show that Schnirelman's property implies hoo = 
on T^. The proof proceeds ab absurdo. 

Let us assume that hoo ^ 0. Using the notation f~{x) = min(/(a;), 0), this entails 
/t2 h^{^)dx = —A < (we recall that h^o is in L^(T^), and Jj2 hoo{x)dx = by construc- 
tion). Then, the theorem of dominated convergence implies that /j.2 h^ j^{x)dx < —A/2 
for m large enough, hence the open set -Em.-A/a = {x E \ hoQ,m{x) < —74/3} has a non- 
zero Lebesgue measure. Now, for any e > 0, lemma 1 implies that for all N > N{m,e), 
X e -E^ _^/3 =^ hN{x) < —A/3 + e, and therefore 

f ^Nh^(.) < £(^^__^/3) e^(-^/3+^) . (28) 

(We choose e < ^4/3 for the sake of the proof). 

On the other hand, the compact set = {a; G | hoo{x) > 0} has non-zero 
Lebesgue measure, since /p^ hoo{x)dx = A. We want to estimate the size of its subsets 
F;v,+ — ^ ^+ \ hN^x) > 0}. The weak-* convergence h^o entails that for N large 



21 



enough (say, > No), A/2 < Jp^hi^{x)dx, itself less than or equal to /p^^ /i7v(a;)(ia;. 
Since all the potentials /ijv are uniformly bounded by M, we obtain the lower bound 
t^{FN,+) > A/2M for all N > No. Therefore, 

/ e^''^(^) dx> [ e^''^(^) dx > C{Fn,+) > A/2M. (29) 
We now combine the equations (28,29) and obtain 

m > max(iV(^, e). No), /'^''fl'^l > e^^^^'-^ (30) 

^E^,-A/z^N{x)dx 2MC{Em,-A/3) 

Such a sequence of Husimi densities {i^Ar} obviously violates Schnirclman's property, since 
the above ratio diverges in the limit N ^ oo, instead of converging towards 
This concludes the proof ab absurdo. 

We have proven that for any sequence {Hj^j} weak-* converging to C, the only possible 
accumulation point (in the weak-* topology) of the corresponding sequence {pn} is C. 
Since this sequence stays in a compact set, we deduce that it converges weak-* to C 

Theorem 1 For any sequence of Husimi densities {H .(n) }NeN weak-* converging to the 

Lebesgue measure C on the torus in the semi- classical limit, the corresponding densities 
of zeros {p,(iv) }N£n also equidistribute in the weak-* sense. 

This theorem concerns Husimi densities defined on the torus, but it can be easily gen- 
eralized to sequences of Husimi densities (and the associated constellations) living on 
compact Riemann surfaces of any genus, using the formalism sketched in section 2.1.1 [6]. 
The theorem has also been recently generalized [25] to holomorphic sections sn{z) on the 
powers L®^ of a positive hermitian line bundle L over a compact Kahler manifold (X, u) 
of arbitrary dimension; the analog of the Husimi density is the hermitian metric of the 
section Hs^{z, z) =^ ||sjv(-2)|P; the zero set of sjy defines a (1, l)-current on the manifold, 
which is shown to converge weak-* to the Kahler form cu in the limit N ^ oo a,s long as 
the Husimi densities become uniform. 

In the case of the 2-sphere, the asymptotic equidistribution of chaotic eigenconstella- 
tions had already been noticed from numerical calculations [22, 23]. The data actually 
showed a more precise phenomenon: not only do the zeros spread throughout the whole 
phase space, but they also seem to repel each other at distances of order 1 / \/iV (see sec- 
tion 5 for a more quantitative statement of this phenomenon). Is such a repulsion another 
consequence of Schnirclman's property? 

4.2 Further considerations on Schnirelman's property 
4.2.1 Fourier coefficients 

In the toral geometry, the weak-* convergence of a sequence of positive normalized mea- 
sures Pn to Poo is equivalent to the convergence of each Fourier coefficient (eq. (21)); we 
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have just proven that Schnirelman's property for a sequence {i^Af} imphes 

e z^ ^ 5k. (31) 

(5fe is the Kronecker symbol as in eq.(23)). We would hke to know better the rate of 
decrease w.r.t. N of these Fourier coefficients, as well as their dependence on k, for N 
fixed. Apparently, this kind of information is not easy to extract from weak-* estimates. 
However, by formally looking at the relation 

Hn{x) = exp X: ^ e^-'^-j , (32) 

(implied by eqs.(22)-(23)), we expect the following: 

Conjecture 1 For any sequence {H (n) } Nen with the Schnirelman property, the k ^ 
Fourier coefficients of the constellations decay as p^r,fc — o{N~^) for N oo. 



We are presently unable to prove such an assertion, even with N ^ replaced by some other 
definite o(l) function. 

4.2.2 Upper bound on the potentials 

Alternatively, we would Uke to control more precisely the semi-classical properties of the 
electric potentials h]sf. Via a Parseval-Plancherel formula, we control their norms: 



/ 



I n l2 

E 



\PN,k\ 



< i: ^+ 1: 



2 



Of^\k\<Ki ^ \k\>Ki ^ 

since all Fourier coefficients arc normalized. Adjusting the cutoff Ki and using the con- 
vexity property < r < s < ||.||s, we prove the following: 

Theorem 2 Let {ijjN{z)}N be a sequence of Bargmann functions s.t. their densities of 
zeros {pn}n weak-* converge to C Then, for any 1 <s <2, the corresponding potentials 
{Hn} tend to zero in the U norm, as N ^ oo. 

In the general case studied by Schiffman and Zelditch [25], the above convergence was 
proven for the norm. 

To obtain pointwise information on the potentials, we can use lemma 1 with the 
knowledge that h^{x) = 0: this yields 

Ve, 3N{e) s.t. > 7V(e),Vx e T^, hN{x) < e. (34) 
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We then naturally tried to estimate the rate of increase of A^(e), as e — > 0. In other words, 
we searched for a universal increasing function f{N), such that the functions f{N)h]\r{x) 
arc uniformly bounded above by a constant independent of A^, as long as {H^} have 
the Schnirelman property. The following example of sequences of Husimi densities shows 
that we cannot do better than f{N) = 1 in general. At the same time, it shows that 
Schnirelman's property does not preclude high degeneracies of the zeros. 

4.2.3 Example of singular Husimi function with the Schnirelman property 

We present below sequences of Husimi densities which look quite singular as far as their 
smoothness is concerned, although they weak-* converge to C These Husimi functions 
are built by rescaling and taking a certain power of the elementary H^(x), so their zero 
constellations will consist of square lattices, with each zero multiply degenerate. For any 
couple of integers {Ni,N2), we consider the Husimi function in Tijv, {N = N1N2) defined 
as 

H^{x) e^^'^^^i-) . (35) 

On the one hand, the 1/Ni periodicity ensures that the Fourier coefficients H^k and 
PN,k vanish unless A^i divides both entries of k. Therefore, if we consider a sequence 
{N — NfN2} s.t. A^i — » 00, Schnirelman's property holds for the functions {H^} (and at 
the same time, Npi^i^k for all k). 

On the other hand, all invariant functionals take the same values for Hn and the 
completely degenerate state of order N2 (cf. fig. 2b, right). If N2 gets large, we can use 
the asymptotic results from the last fine of table 1, replacing N by N2. In that case, 
the Husimi densities Hj^ asymptotically look like delta peaks concentrated on the lattices 
dual to their constellations, which explains the large deviations of the invariants from the 
totally flat case (i.e. H = 1). 

The electric potentials of the constellations are given by h^lx) = ■^h{Nix)^ bounded 

above by M/Nf. If N2 increases much faster that Ni (e.g. N2 ~ e^O' then a function 
f{N) as defined above has to increase slower than (e.g. ~ log^ N). There is therefore 
no universal strictly increasing f{N). 

At the same time, we notice that such a singular behavior of the invariants, and of the 
shape of the densities, comes along with a high degeneracy of the zeros. In that sense, the 
strict repulsion between the zeros of chaotic eigenconstellations, conjectured in [22, 23], 
cannot be deduced from Schnirelman's property. 

4.2.4 Example of non-ergodic Husimi functions with equidistributed zeros 

We exhibit a sequence of Husimi functions showing that the converse of theorem 1 is 
false. We use constellations in the shape of deformed lattices, defined as follows. We 
start from a smooth separable density on T^, p{x) = Pq{q) Pp{p) , with the normalizations 
lo Pi{v)dv = 1, for i = q,p- Separability allows us to integrate this density, i.e. to change 
coordinates x 1— >• M{x), s.t. 'p{x) = det-^(a;) (we just take Mi{v) — Jq pi{v')dv' for 
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Now, we build the Husimi function Hp^N^x) with zeros at the N = points Xij = 
A^-^ (i/M, j/M) , through 

M 

h-p,^{x) logH-Mq,p) E h {q-N-\i/M),p-N;\3/M)) . (36) 

We want an approximate formula for H-pj^{x), in the limit M oo. For that, we 
transform the right-hand side by the Poisson summation formula at fixed x, which yields 
two types of terms. The zero Fourier coefficient (Weyl's term) is M'^h*p{x), whereas the 
others (oscillatory terms) read 

/ h{y)p{x + y) e-^i^rMfe-AAC^z+x) ^2^^ Q^kei?. (37) 

For /c 7^ fixed and M +00, the above integral is dominated by the contribution 
near the singular point j/ = 0. Expanding M around x and summing over all k ^ 0, we 
(formally) obtain 

hp,N{x)-M^h*p{x) ~ logH^(r(x)){MAfg{x)+r{x)MAfp{x)), with r{x) = (38) 

{H-^ij^ is defined in equation (17)). The convergence of the Poisson series to this function 
is uniform as long as one stays away from the Xi j (notice that their positions change with 
M !). Anyway, since we subsequently exponentiate this function, we need not pay too 
much attention to problems of convergence near the zeros. In the large-M limit, we thus 
obtain a "WKB-likc" expression for H-p^]^{x): 

H-p,n{x) ~ exp(iV/i * -pix)) H^ir{x)) {MJ\f^{x) + T{x)MMp{x)) . (39) 

Locally, this oscillatory function vanishes on a rectangular lattice of spacing oc 1/M, since 
the parameter t{x) varies slowly compared to this spacing: the large-M constellations thus 
resemble deformed lattices. 

We now let the smooth density p itself depend on N: we take the sequence A/jv(a;) = 
{q+XqSm{2nq),p+XpSm{2TTp)), which yields p^(a;) = (l-|-27rAgCOs(27rg'))(l-|-27rApCOs(27rp)), 
where the coefficients depend on A^. If Ag, Xp vanish in the limit N — > 00, then J)^ 1 
and t{x) i uniformly on so that the oscillatory factor has asymptotically the local 
average 1/GM[H^]. 

Meanwhile, the convolution under the exponential, h * p{x), yields the trigonometric 
polynomial Pn(x) — 2[(7rAqCos(27rg)-|-l)(7rApCos(27rp)-|-l) — 1]. Then, if the decrease of 
either Aj is slower than l/N, the factor exp{NPiy{x)) is singular in the semi-classical limit. 
For instance, if both Aj decrease slower than it gives a peak (an artificial point scar) 
at the point (1/2, 1/2); whereas in the marginal case Aj = Ai/N, with Aj non-vanishing 
constants, that factor is asymptotically A^- independent; Hp^,N then converges weak-* to 
exp(— 27rAg cos(27rg) — 2nAp cos{2np)) / GM[H^], a non-uniform density. 

So, Schnirelman's property does not hold for such sequences Hp^^N, although pjv ^* C 
However, this counterexample is just consistent with the converse of conjecture 1. 

In figure 4 we display Husimi densities of the type described above, for various values 
of A" and Aj, together with the corresponding invariants. 



25 



4.3 Phase of the Bargmann eigenfunctions 



From the semi-classical equidistribution of the eigenconstellations, one can get some esti- 
mates about the phase of the Bargmann wave-function, which appears as a complementary 
information to the Husimi density (see eq.(5)). For this purpose, we compute the deriva- 
tive of \ogH{x) w.r.t. the holomorphic variable z, and get the following representations: 



27rNz = 



d_ 
dz 



log H{z,z) 



= N[{d,h)*p]{z,z) 
= -iNV2 J2 



Pk 



-^2mk.x 



(40) 

(41) 
(42) 



We will use the second and third equations to prove the following theorem. 



Theorem 3 Let {ipNiz)}^ be a sequence of Bargmann functions s.t. their densities of 
zeros {pn}n weak-* converge to C Then, for any 1 < s < 2, the functions 2^^(z) 
tend to the function z in the norm, as N ^ oo. 

To prove the theorem, our strategy is to transform the series (42) into a finite sum, up 
to a small remainder, and then use the convergence of the individual Fourier coefficients 
PN,k (cf. eq.(31)), as in the proof of theorem 2. However, due to the stronger singularities 
of the logarithmic derivatives, we now need to first regularize the functions dzh * pn by 
Gaussian convolutions dxix) = K"^ e~'^^ ^ . In the limit K — > oo, 5k converges to the 
identity kernel S. For finite K, the deviation from the identity is given by the series 



[dzh - dzh * Sk] {q, p) = ^2 ^ 



n,m& 



Q-nK^[{q-mf+(j>-n)^] 

{q — m)+ i{p — n) ' 



(g,p)^z2. (43) 



Using the normalization of pat and the equation above, we obtain the following estimates, 
valid for any fixed 1 < s < 2: 



T2 



[pN * {dzh - dzh * 5k)\{x) 



dx < 



f2 



< {2ns) 



[dzh-dzh*SK]{x) 
r(l-s/2) 



dx 



s/2'_ 



(44) 
(45) 



The first inequality is just due to convexity (it is equivalent to the triangular inequality 
for the norm); the right hand-side can then be calculated exactly, leading to (45). The 
divergence as s — > 2 is obviously due to the poles of V'VV'- Once we select the power s, we 
can adjust the cutoff parameter K large enough to make ||pAr * {dzh — dzh * Sk)\\s small. 
On the other hand, the smoothed term can be now estimated: 



\pn * dzh* Sk\\s < Wpn * dzh* Sk\\2 



(46) 
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< 12: T^K^re— ^1 (47) 

In the last inequality, we introduced a cutoff Ki, and estimated the remainder of the 
series J2\k\>Ki using the uniform bound \pN,k\ < 1- For Ki/K large enough, this remain- 
der is small; finally, the finite sum J2\k\<Ki vanishes in the limit N ^ oo due to the 
equidistribution of zeros. 

Notice that the theorem extends to values of s in the interval (0, 1), due to the convexity 
property < r < s < oo =^ ||.||r < | Mis- 
Interpretation 

In view of the formal analog ih{ilj' / ip) (q) p — the momentum — in the 1-d 
Schrodinger representation, we can define tt{z) =^ ^^{^) ^ the quantum local momentum 
of the Bargmann function (whereas z is the classical symplectic conjugate variable of z). 

The Cauchy-Riemann relations link the estimates 7r{z) ~ z of theorem 3 to the vari- 
ations of both the modulus and the phase of the Bargmann functions: 



l=(tj-{z)-2nNz^ = {d,\og\ij{z)\-7rNq) + i{dp\ogmz)\-7rNp) (49) 

= -{dpeirgijj{z)+TTNq)+i{dqaxgilj{z) -nNp). (50) 

These two equations correspond to interpreting the zeros respectively as 

• either point-like electric charges, generating a potential —log\ip{z)\, balanced by a 
uniform charge distribution of potential ttNx^ /2 (the total electric charge on the 
torus vanishes). 

• or vortices of magnetic flux: the corresponding vector potential is the gradient 

— * 

Varg'^, so that each vortex carries a unit of flux = — 27r. The additional vec- 
tor potential 'kN corresponds to a uniform magnetic field B — 2'kN (the 
total magnetic flux on the torus vanishes). 

It seems difficult to assert the approximation 'n:{z) ~ ^ in stronger topologies than 
the L" estimates of theorem 3. Generically, the strict equality ti{z) = z can only hold 
at isolated points on the torus, due to the analyticity of ip{z), or equivalently due to its 
phase. Indeed, the right-hand sides of equations (49,50) have quite different behaviors if 
we try to set them equal to zero. The equation on the modulus has the obvious smooth 
solution log |-?/'(2;)| = ttNzz, which corresponds to a strictly uniform Husimi density, i.e. 
the formal Schnirelman limit. By contrast, the equation on the phase reads 



Varg^(.)= (51) 
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which has no solution, since the corresponding magnetic field V x ^nNq j ~ "^ttA/" 

does not vanish. The singular behavior of the phase of ip{z) in the semi-classical limit 
can be therefore interpreted as a 'struggle' to achieve as well as possible the equality (51); 
such a singular behavior implies a distribution of phase dislocations (i.e. of zeros of ijj) 
all over in the semi-classical limit. 

At first sight, the situation looks quite different in the integrable case, say for eigen- 
states of a time-independent Hamiltonian which we express in the complex torus coordi- 
nates as Ti.{z,z). Assuming that TC is an analytic function and that the eigenvalue E is 
a regular energy value, the WKB form (1) holds for semiclassical eigenf unctions in the 
Bargmann representation [21]. Then, in the h ^ limit, h{logiJjy{z) ~ S'{z) solves 
'H{z,S'{z)) = E, the Hamilton-Jacobi equation in the z variable: hence S'{z) = yE{z), 
where {y — yE{z)} is the classical complex energy curve Eg = {l-i{z,y) — E} solved for 
the conjugate momentum of z. These asymptotic forms hold outside the anti-Stokes 
lines, i.e., the curves where any two real parts of the (multisheeted) function ysiz) match. 
All in all, the dynamical result for an integrable eigenfunction is then 

7r{z) ~ yE{z), outside anti-Stokes fines, for ^ — > 0. (52) 

Precisely, the quantum local momentum of the Bargmann function 7r{z) tends to some 
branch of the classical local momentum function, ysiz) in the ^-representation. 

However, the perspective changes if we seek the asymptotic behaviors specifically 
obeyed in the classically allowed part of phase space (which is the only scmiclassically 
meaningful region, where the Husimi density will not become negligible). Here this region 
is the real energy curve Eg = {ysiz) — z} (locally), on which eq.(52) implies 

^{z) ~ ^ in the classically allowed phase space, for — 0. (53) 

We now compare with theorem 3 for an ergodic situation: then the classically allowed 
region is the whole phase space, so we realize that the same asymptotic result (53) has 
become generalized from integrable to chaotic situations, in which case it moreover applies 
almost everywhere in phase space — albeit in a weaker sense (L*) than before (pointwise). 
The contrast between implementations of eq.(53) for integrable vs chaotic dynamics is 
illustrated by color plots for 7r{z) (fig.5 top); we add that the same overall appearance is 
universally shown by all such plots made for equidistributed zeros (disregarding the precise 
distribution of the singular points themselves), whereas it is different and case- dependent 
for zeros concentrating on curves. 

Upon a purely formal integration, eq.(53) yields 

" logiP{z) r-. ^j' z'dz' (54) 

but this is as inconsistent as eq.(51) for specifying Q\ogtlj{z), since the 1-form Q{z'dz') is 
not closed. Hence the phase of ip{z) will stay undetermined pointwise; however, eq.(53) 
specifies as universal its variation pattern when zeros are equidistributed, and this is 
displayed by color plots for ip{z) itself in fig.5 bottom. 
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4.4 Stronger estimates for lattice eigenstates 

We can strengthen the above estimates for a particular class of eigenstates of quantum 
cat maps [38], namely eigenstates for which the constellations form lattices on the torus 
(such a lattice must be invariant under the classical cat map) . The construction of these 
eigenstates is performed in [12, 42]. For a given classical cat map S, such states exist only 
for some particular values of the inverse Planck's constant A^, in which case they are quite 
scarce (they generate a subspace of T^at of dimension small compared to A^). Moreover, 
their quasi-energies are usually degenerate. On the other hand, the regularity of the 
constellations allows us to estimate more precisely the invariant functionals, as well as 
the phase variations. Indeed, an alternative characterization of such a lattice state is that 
its Bargmann function can be written as a single Jacobi thcta function [42]. For instance, 
if the eigenconstellation is the lattice on generated by the two complex numbers [^1,^2] 
ordered s.t. '^{r =^ V2/V1) > 0, then the corresponding Bargmann and Husimi functions 
read 

where the elementary functions x(2;|r), -f^^M defined in equations (17). ^^{z) is then 
quasiperiodic (and is periodic) w.r.t. the lattice [^1,^2], which is a stronger property 
than eq.(2) (the Husimi density of such an eigenstate is displayed in fig. la, top right). 
Precisely, for any point v — nvi + mv2 of the lattice, we have 

tP{z + v) = (_i)("+"^+«"0e-^M'+2'^^^^^(^) (56) 
^^(z^v) = ■^iz)^2T:Nv. (57) 

A crucial property of these lattice eigenstates is that for a given cat map S, the values of 
the modulus r stay inside a compact domain of the upper half-plane, namely a rectangle 
|3'^(t)| < 1/2, 1/2 < S>(t) < C5, with Cs a constant independent of (see Appendix B). 

The norms and related functionals introduced in section 3.2 are invariant through 
a rescaling of the variable z, so they take the same values for H^{x) of equation (55) 
and H^(^T--j{x). Since these functionals are smooth functions of r (cf. App. A), they are 
bounded uniformly w.r.t. N for all lattice eigenstates ip of the quantum operators Us 
associated to a given cat map 5" (the bounds depend on 5") . 

Similarly, for such lattice states we get fine estimates of the logarithmic derivative of 
ip{z), using the quasiperiodicity (57). We obtain from eq.(55). 



^(z) - 27rNzf£z = \^idg + idp)\ogH^^,^{q,p)[dqdp. (58) 



As r stays inside a compact domain, the integral of |9^1ogif^(T-)|^ on T^- is bounded. On 
the other hand, we show in Appendix B that l^il"^ = 0{\/N), so we may strengthen 
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theorem 3: 



z) -z\\s< with Cs,s indep. of N. (59) 



On the Fourier side, to a given cat map 5" is associated a constant Ks s.t. for any 
lattice eigenstate ipN, the Fourier coefficients of the constellation (which are supported by 
the dual lattice) have the following property [42] : 

VA; 7^ 0, \k\ < Ks^ =^ PN,k = 0. (60) 

Notice that the above equation trivially implies that conjecture 1 holds for a sequence of 
lattice eigenstates, since pjvA; = for large enough. 



5 Statistical model 

In section 3.3, we observed that the values of the different functionals for eigenstates 
of quantum chaotic maps were surprisingly close to their average values over a certain 
ensemble of random vectors oi Hn, which indicates that a chaotic eigenstate 'looks like' 
a random state, at least in a certain sense. This remark is linked to the various ran- 
dom matrix conjectures in quantum chaos [32], since the eigenstate of a random matrix 
is a random state. We will not try to justify these conjectures in the following, but 
rather describe further the relevant random states and their properties, especially in the 
Bargmann-Husimi-stellar framework, in order to compare to them the corresponding 
quantities computed for eigenstates of quantum chaotic maps. 

Such statistical models were already studied in alternative phase spaces, namely the 
2-sphere and the plane [27, 29, 23, 26, 28], where the Bargmann functions arc respec- 
tively polynomials or entire functions in z of controlled growth. The main observations 
concerning the statistical constellations were the following: 

• first, the zeros are on average equidistributed over the phase space, or the classically 
allowed part thereof. This property matches the equidistribution we proved for the 
eigenstates with the Schnirelman property (theorem 1). 

• second, zeros at short distance (approximately \/h) tend to repel each other [23, 26, 
28]. Note that this typical distance coupled with the finite phase space volume imply 
a certain rigidity of the constellation. This phenomenon seems present also for eigen- 
states of quantum chaotic maps; nonetheless, we already argued that such a repulsion 
cannot be explained by Schnirelman's property (cf. section 4.2.3). To characterize 
precisely the rigidity of the random constellations, one could study the statistics of 
their Voronoi tessellations (such tessellations are shown in figures la and lb, bottom 
rows, for chaotic eigenstates). This approach was used for instance to characterize 
eigenvalues of large complex random matrices [33]. We did not investigate in this 
direction in the present article. 
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Symmetries 

The several models we will introduce correspond to vectors with different symme- 
try properties; symmetries play a fundamental role in the spectral properties of random 
matrices [32], and their relevance in the description of random states has already been 
noticed [26, 28]. 

The classical maps we study are all invariant under parity {q,p) i— > {—q,—p), and 
the quantizations we consider then yield either even or odd eigenstates; by linearity, the 
Bargmann functions of these eigenstates will also be even or odd functions of z. As a 
consequence, we will consider models of even random states (the treatment of odd states 
is very similar, and we skip it). 

Besides, the classical maps can also have anti-canonical symmetries. The baker's 
map is invariant under the reflection {q,p) ^ {p,q), and the simplest 'cat' map usually 

/ 2 1 \ 

considered (i.e. the matrix ^ ~ { ^ 2 V time-reversal symmetry {q,p) 

(g, —p). On the quantum side, these symmetries manifest themselves on the coefficients 
of eigenvectors. For instance, the time-reversal symmetry corresponds to eigenvectors 
with real Schrodinger coefficients, giving real Bargmann functions (i.e. iIj{z) = ipiz)). 
Such a property has to be incorporated in the statistical model as well [26, 28]. 

It is possible to build chaotic maps on the torus with no anti-canonical symmetries 
(some generalized baker's maps, or for instance the cat map S' of eq. (13)). However, 
these systems are in practice much more hyperbolic than the simpler ones cited above, 
so that their semiclassical properties (e.g. scars) may appear only for very large values 
of N, which is inconvenient numerically. We nevertheless start our study by considering 
an ensemble of random vectors with no symmetry whatsoever; this model yields tractable 
analytical results, not all of which can be generalized to more complicated ensembles. 

5.1 Model with no symmetry 

We define a random state in Hn,o — ^ without loss of generahty), as: 

Af-l 

V'{a} = Yl %kj)jV,0, (61) 
j=0 

where the Schrodinger coefficients aj are random independent Gaussian complex variables 
with (%)^ = 0, (ajakj^ = Sjk/N. The states \qj)N,o-, j = 0,...,A^ — 1, form the 
orthonormal basis of position eigenstates of Hn^ [42] (but the statistical ensemble is 
invariant through any unitary change of basis) . 

The vectors ip^a} are not normalized a priori, but their square norm 712 — ||'0{a}||^ has 
the distribution law 

increasingly peaked around ^2 = 1 in the limit N 00. With this fact in mind, we can 
still calculate the typical values of the different functionals introduced in section 3.2. 
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The distribution law of the Bargmann function at a given point Zo is the Gaussian 



— e (fib with width 

Y N-1 



a 



Zo 



N j 



E l(^ok,)^,or. (63) 



i=o 



This width is actually a particular case of the correlation function (^^jJ{zl)^lJ{z2)j^ which 
happens to be the crucial quantity of the statistical model. We give its value for an even 
integer (the formula for odd is slightly more complicated, but has the same large-A^ 
behavior) : 

AT-l 

(V'(^i)V(^2))^ = N-'J2i^i\Qj)N,o{z2\qj)j,,o (64) 

j=0 

= N'-^ N,o{zi\z2) N,Q 

^ ^2nm,z, g^^^]y^l±Jl I i^/2) OsiinN^^^ \ iN/2). 

y2 V 2 

The above formula for overlaps of coherent states on the torus is compatible with the 
corresponding estimates derived in [14, §8.2] and [15, §1.4.2]. We deduce from this that 
the distribution law of the Husimi function at a given point x = {q, p) is the exponential 

1_ Q-H/a. ^ e^inNq I iA^/2) 63(7: Np \ iN/2), (65) 

and the limit dx ^ 1 holds exponentially uniformly on as — > 00. Notice that 5"^ is 
also the square-norm of the torus coherent state |^)Ar,o, up to the factor N e^'^^^^. Prom 
this pointwise distribution, we can derive the following averages, leading to the typical 
values of several invariants as given in table 1, 

(J^^H{xfdx)^ = 2^3(0 I iiV)2 ~ 2 + 8e-'^^ 

/ H{xydx[ Hiyfdy) ~ r{l + s)r{l + t) \l + ^ T - [^F.i-s, -t;l;u) - 1] du] 

(I logH{x)dx) = --f + 2logri{iN/2)+TiN/12 7-2e-'^^ 

f H{x) log H{x) dx) = 1 - 7 + 2 Osinq \ iN/2) log Osinq \ iN/2) ~ 1 - 7 + 2 € 
' n Jo 

{•zFiia.b; c; z) is the hypergeometric function [2, vol. 1, chap. 2]; the function rj^r) is 
Dedekind's modular form defined in appendix A; 7 is Euler's constant). The second 
identity allows to derive the typical deviations from average values for the L^-norm, the 
geometric mean and the entropy. 

The average sup-norm of a random state cannot be calculated by the same techniques 
as the other invariants. We can nevertheless obtain an upper bound for it (see Appendix 
C), as given in table 1: 

' |i/||oo)^<21ogAr. 
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This result (which is the phase-space counterpart of similar estimates for sup-norms in 
the Schrodinger representation [34]) actually concerns a constrained ensemble of random 
vectors, namely the set of complex vectors on the unit sphere <S^ of H-Nfii equipped with 
its standard measure: 

^ N-l 

p{ao, aN-i) = Vol(^iV) ^(^ ~ ^ I'^j'^)- 

According to equation (62), in the scmiclassical regime this (microcanonical) ensemble 
must yield the same average values as the (canonical) Gaussian ensemble studied so far, 
so we expect the above upper bound to hold for the Gaussian case as well. 



5.1.1 Statistics of the constellations 

Alternatively, it is possible to extract statistical properties for the constellations of the 
random Bargmann (or Husimi) functions, as was done in [26, 27] for a spherical phase 
space. The methods used therein can be transposed to the torus. Indeed, Hannay [27] 
shows that any joint probability function pk{zi, Z2, ■ ■ ■ , Zk)(Pzi . . . (Pzk of the zeros of ran- 
dom Bargmann functions can be computed from a unique 2k x 2k correlation matrix, 
whose entries are the functions (ip{zi)ip{zj)'j ^ and their derivatives w.r.to Zi or zj. Notice 
that both the Gaussian and microcanonical ensembles of random vectors exactly yield the 
same statistics for the zeros, which are independent of the normalization. 

1-point function 

The first relevant quantity is the average density of zeros on the torus, defined as 
Pi{Z) = (^N~^ J^fJo' ^{z — ^j)}j^-> where the Zj are the zeros (ordered arbitrarily) of a 
sample ip^a}{z). The S functions are two-dimensional, and periodicized w.r. to T^. These 
notations will apply as well to the pair correlation function defined below. Hannay's 
formula reads 

p,{z)d'z = (69) 

^Pi{q,p) = I + ttA^E - L at ,^A cos{27rNmq)+cosi27rNmp)] (70) 

^1 smh(7rArm/2) 

■kN 

~ 1 - h( jv/2) ^^"^^^^'^^^ cos(27rA^p)] as ^ oo. 

Therefore, in the scmiclassical limit, the zeros of random Bargmann functions are equidis- 
tributed over Tc- The deviations from perfect equidistribution {pi{x) = 1) are uniformly 
exponentially small, and periodic w.r. to the square lattice of side 1/N: this seems re- 
lated to the sum rule (7), which prevents the distribution of the constellations from being 
completely translation- invariant . 

Derivative of ^ at the zeros 

The Husimi density is very oscillatory in the scmiclassical limit, due to the dense dis- 
tribution of its zeros. In the vicinity of a zero Zi, it can be approximated by H^{z, z) ~ 



33 



\{z — Zi)il)' {zi)\^ e^'^''^''^'''. Therefore, the quantity f{zi) ^= \il)'{z-^\^ c^'^'^^''^''^ seems a reh- 
able scalar to measure the local strength of if^ around Zi as built by all the other zeros, 
and it fluctuates much less than H^{z, z). 

In our statistical framework, using the joint probability of the Bargmann 

function and its derivative at a point Zi, we can derive the distribution law of A^~^ J2iLi fi^i) 
for {zi} the zeros of ip: in the semiclassical limit, this law happens to be the so-called 
"1/2-Poisson law" 

"df. (71) 



I. f.-f/-2nN ^ 



(27rNy 
2-point function 

The exact pair correlation function P2{zi, z^) — (^N"'^ ^ij^j ^(^i~-^j)^(^2— -^j))^ would 
already be too lengthy to write down, even in this model without symmetry. However, 
the formulae get much simpler in the semiclassical limit, i.e. by neglecting exponentially 
small corrections: we then recover an expression very similar to the one applying to 
the spherical phase space. Actually, we use the fact that the fundamental correlation 



function {^^[zi)ip{z2)j ^ is given semiclassically by the kernel e^^^^^^'^, uniformly when 
the separation 5z — z\ — Z2 stays inside a compact square < 1/2 — e, \Sp\ < 1/2 — e (in 
such a square, the thcta functions in eq.(64) converge uniformly to 1). In this regime, the 
simplicity of this kernel yields for p2 a formula both translation-invariant and isotropic, 
identical to the one found for the sphere [27, 26]: 



P2(xi,X2) ~ g{^7rN/2 \Sx\) as ^ oo, (72) 
, , def (sinh^ -I- r^) cosh — 2r^ sinh 

= ^h^hV ■ ^ ' 

The function g{r) is displayed on figure 6 (left). It starts quadratically near the origin, 
and converges to 1 ~ Pi{xi)pi{x2) as soon as l^i — ^ 1/y/N. The zeros (interpreted 
as interacting particles), repel each other at short distance, and become uncorrelated at 
distances larger than the mean spacing 1 / 

In the case where \Sq\ or \dp\ is near 1/2, we cannot use the kernel e^'^^^i^^ g^j^y more, 
but a direct estimation of the different terms contributing to p2 leads to the large-iV 
uniform value ^2(2^1, 3^2) ~ 1, as was expected from the shape of g{r). The periodicity of 
P2 w.r. to then yields a uniform semiclassical approximation for it valid everywhere. 

Fourier coefficients 

In a series of articles [23. 24, 26], formula (72) was compared with numerical compu- 
tations using eigenstates of quantum chaotic maps on the sphere or on the torus, and an 
excellent agreement was found. However, averaging over many eigenstates was necessary 
to recover the shape of the statistical result. Moreover, such a direct comparison yields 
information about the short-distance correlations between the zeros, whereas we are also 
interested in the global properties of the constellations. In fact, the semiclassical proper- 
ties of eigenf unctions are certainly not given by the precise position of a given zero, but 
rather by interferences between a large number of them. Consequently, as we already 
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explained in section 3.5, we think that some relevant semiclassical information may lie 
in the Fourier coefficients of the density of zeros rather than in their individual positions 
(see section 6 below). We can actually estimate these Fourier coefficients in our statistical 
framework. 

The formula for the average density pi{x) shows that (pfc)^ = up to exponentially 
small corrections. 

More information is contained in the second moments of the coefficients. We are thus 
led to the form factor of the random constellations, defined for any A; e by 



\Pk 



'n N 



\-+ I d'x, d'x2 e-2-^-(-i--^) P2{x,,X2), (74) 

where p2 is the pair correlation function corresponding to random states in T-Cn,o- Using the 
asymptotic formula (72) for p2 and integrating over the angular variable, then integrating 
by parts, we obtain the following integral and series representations, valid as iV — > oo up 
to exponentially small corrections. To this order, the isotropy of p2 in (72) implies the 
isotropy of the form factor: 

{\Pk\')^ ~ 6, + N-'F2{\k\/^/N) (75) 
with F2{k) = -27rK^ y J2{n^STry) (1 - cothy) dy (76) 

n>0 ^■ 

where J2 is the Bessel function and C,{n) — J2j>oj~"' the Riemann zeta function. Rigor- 
ously, ^Ipfcp^^ is only defined for integer k, but the limit N ^ 00 involves the values of 
F2{k) for K G M+. We were not able to find a simpler expression for the rescaled function 
F2, which is shown on figure 6 (right). The curve F2{k) looks vaguely similar to the 
curve g{r) of figure 6 (left), but the two functions have different behaviors near the origin. 
Whereas g{r) ~ for small r, the series (77) shows that for fixed k ^ 0, 

^J{\Pk\^)^ ~ 7ryC(3) \kf N-^l^ when TV ^ 00. (78) 

Therefore, the decay property stated in conjecture 1 is amply fulfilled on average for 
sequences of random states. 

Notice that we also obtained here the typical dependence of as a function of k for 
fixed iV, of which we had no idea from our deterministic point of view of section 4. 

The linear relation (23) in Fourier space between the density of zeros and the logarithm 
of the Husimi density provides us with some statistical information about the latter: 

|/i,|2\ ^AT-sM^ as iV^oo. (79) 

This scahng function is also plotted on figure 6 right (dot-dashed curve). In particular, if 
one fixes 7^ A; e Z^, then 

|2\ 



A^' (|/ifcr)^ ^ C(3) as N^oo. (80) 
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5.2 Even-parity model 

The statistical model we have just studied had the advantage of translation invariance 
and isotropy, up to exponentially small terms in the semi classical limit. However, the 
chaotic maps on the torus that we studied numerically have the parity symmetry, and their 
eigenstates are either even or odd, in both the Schrodinger and Bargmann representations. 
This property then has to be incorporated in the statistical model. We will only write 
down the statistical properties of random even states, and to simplify the formulas we 
restrict ourselves to the case even. 

The correct ensemble is built in the following way: we take a random Bargmann 
function ip^a}{z) from the ensemble with no symmetry (61), and consider the symmetrized 
function 

^even^^^) d^f ^{a}{z) + lP{ay{-z) ^^^^ 



N/2-1 



ao(2;|go)iv,o + — — j^^{{z\qj)N,o + {z\qN-j)N,o) + \/^aN/2{z\qN/2)N,o 
.7=1 v2 



N/2-1 



= ho{z\qQ)Nfl + hN/2{z\qN/2)Nfi+ hj({z\qj)N,o + {z\qN-j)N,o), (82) 



so the new random independent Gaussian complex variables 6j have the variances (l^op)^^ 



hN/2?)^ = )^ = for j = 1, . . . , N/2 - 1. For this ensemble too, all in- 

teresting quantities can be evaluated by use of the correlation function 

(^(^0^)7'^ = (^(^0^)^ + {i^{z,)J{^) ^. (83) 



even 



Therefore, we have the uniform semiclassical approximation (jjj(zi)ijj{z2)j ^ ~ 2 cosh (27rA^2;i 2^2) 
as long as both {xi — X2) and {xi + X2) stay inside a square of side 1 — e centered on the 
origin (in this square, this formula holds up to exponentially small corrections). 

The distribution laws of the Bargmann and Husimi functions are slightly changed from 
the case with no symmetry. The modification concerns the widths of the Gaussian (resp. 
exponential) distributions: 



= {^l^{z)^{z)) ^ (84) 
a^J^^ = + e-''^^(«'+f') ^3(i7rA^g | \N/2)e^(;mNp \ \N/2). (85) 



The second term imphes that the random even Husimi function has a 'bump' of height 2 
and width ~ 1/VN at the four P-invariant points (0, 0), (1/2, 0), (0, 1/2), (1/2, 1/2) (i.e. 
the points where x — —x modulo T^); at a distance ^ l/\fN from these special points, 
this term is small, and c^™" ~ 1. A consequence of these 'bumps' is that the average 

/ \ even 

square-norm of is now A^-dcpcndcnt: (HV^IP) = 1 + 2/N. The average values of the 



/n 

other functionals similarly differ by terms of order 1/N from their values without parity. 
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1-point function 

The average density of zeros is modified by parity in a manner dual to the 'bumps' of 

the Husimi density, i.e. it decreases in the vicinity of the symmetry points. Indeed, when 

|2\ ^, o^..c,■K/'o-,r/\^l^|2^ 



\q\ < 1/4— e, IpI < 1/4— e and A?" is large, the approximation ~ 2 cosh(27rA'"|2;| ) 

in equation (69) yields: 

pr"(.)~ta„h(.7V.^, + ^;j|^ (86) 

(we show p'^^^lx) for small x in figure 7 left, solid curve). The form of pf™"(a;) in the 
remaining part of can be obtained through its periodicity: pf^^{q,p) = pl^^^{q + 
1/2, p) = Pi™'^(q',p + 1/2). The decrease of the density near the symmetry points can be 
explained as the repulsion between the zeros z and —z which come close to each other 
when they approach the symmetry points. 

2-point function 

The pair correlation function, now defined as p2™°(-2ii -^2) = J2zi=^±Zj ^{^1 ~ 

Zi)S{z2 — Zj)'^ ^, can be calculated by the same techniques as in [27], using the asymp- 



totic correlation function yp(zi)ip{z2)^ ^ ~ 2 cosh(27rA^2;i2;2). The exact asymptotic 
result cannot be written down in a concise way (it depends explicitly on both Zi and Z2). 
However, if xi, X2 are not in the vicinity of the same symmetry point Xgym (i-e. we do 
not have simultaneously \xi — Xsym\ = 0(1/ \/N) and \x2 — Xsy^\ = 0(l/\/iV)), the pair 
correlation function is approximatively 

pr\x,,X2) ~ pr\Xl)Pr'\x2)p2{x,,X2)p2{Xi, -X2) (87) 

where p2 is the 2-point function for random states with no parity. 
Fourier coefficients 

Due to the inhomogeneity of pi, the averages of the Fourier coefficients are no longer 
exponentially small in A^, as was the case with no symmetry. The periodicity of pf"^'^ 
entails that (^Pk)N^'^ = unless kq and kp are both even. In that case, we have 

\ even 1 , , , , / , 

Pk)^ = S, + -F,{\k\/VN) 

roo 

with Fi{k) — 4:K ^/nyJl{2K^/7^y) {ta,nh{y) — 1) dy 

J 

/ \ even 

This provides us with the following asymptotic behaviors of \Pk) (see fig. 8, dashed 
curve) : 

= ^A^- Y^ + C'(A:V^'') for small k/^, (89) 

\ even , i > rC , > 

p,)^ = 0{N-') as ^oc. (90) 
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The form factor is now given by 



A simple Riemann-Lebesgue argument for the second and last terms in this formula 

/ \ even 

implies that that (iPfcP)^ converges to for large {k/\fN), as was the case in the 
former section. On the other hand, using the approximation (87) away from the symmetry 
points, we obtain the following expression in terms of k/y/N: 

Pk?)l^^ = 4 + ^i^2(|A;|/v^) + ^Feo„(^/v^) + 0{l/N% (92) 

where F2{k) is the function given in equation (77). -Fcorr is not known analytically; it is 
bounded and its Taylor series starts by a term of the form ak'^/N with a >0, but we do 
not know if higher terms are isotropic. 

Nevertheless, equation (92) shows that the dominant shape of the form factor as a 
function of k/sfN is unchanged from the model with no parity: we can therefore compare 
this analytical expression to the data for the cat map S' (see section 6 and fig. 9a). 

/ / \ even\ 2 / 1 1 9 \ ^"^^"^ 

Besides, for this model we see that {(pk) ) is of higher order in N~ than {\pk\ ) , 

I \ even 

as functions of so that \Pk)j^ might be difficult to detect (88) in the numerical data 
for large A^. 

On the other hand, if one fixes k and studies the iV-dependence of the coefficient 
(resp. hk), the large- A*" asymptotics are slightly different from the asymmetric model (eq. 
(78) and (80)): 

^^l'>r - -^,{^\m' + ak% (93) 



even 



7V3 

1 / W„N 



We therefore recover the N'^l"^ behavior, but with different prefactors (unless a vanishes, 
of course) . 



5.3 Real even-parity model 



In order to fit the numerical data for the simple cat map S = 




or the baker's 



map, we study a statistical ensemble of real even states (once more, we restrict ourselves 
to even values of A^). This ensemble is still defined by the formula (82), but the bj are now 
random Gaussian independent real variables, with the same variances as in the complex 
case. Models of random real polynomials (without parity symmetry) were studied in 
detail in the literature [29, 26, 28, 30]. 

Due to reality, the distribution laws of the Bargmann function or its derivative do not 
only depend on the modulus l^^l, but also on its phase. However, these distributions are 
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still Gaussian in the variables ^ip, '^ip, "Rip' , '^ip' [28]. On the other hand, the distributions 
of the Husimi function and the averages of its invariant functionals are unchanged from 
the ensemble of complex even random states (former section). 

By contrast, the reality of the hj modifies drastically the statistics of the zeros. 
Since 'ip{z) is real and even, its zeros come either in quadruplets of complex numbers 
{z, —z, z, —z}, or in couples {z^, —z^} situated on one of the four symmetry axes {'^z = 0}, 
{'^z = 1/2^2}, {^z = 0}, {^z = 1/2^2}. The statistics of these two types of zeros are 
quite different, but both are still obtained from the correlation functions (ijj{zi)i/j{z2)'j ^ 



2l 



and {^l){zi)'il){z2)j ^ = {^l){zi)ilj{z2)j ^. Since the hj have the same variances as in the 
complex case, these functions are still given by equations (83,64). 

We use the formahsm of [28] to derive the average density p^™P^^ of complex zeros: for 
\q\ < 1/4 — e, \p\ < 1/4 — e, we obtain semiclassically 

pT^^{x) ~ Gi(v^g,v^p), where (94) 
Gi{Q,P) = ^ /^"^oN./9 {di{Q,P) sinh(g^ + P^) + 2sinh(g^ - P^)[g^ sinh 2P^ - P^sinh2g 

-2gPsinh(g^ + P^) sin(4QP)} 
withcii(g,P) = sinh2g2sinh2p2 + sin2(2QP). 

Figure 7 (right) shows a contour plot of Gi{Q,P) near the origin, in the first quadrant. 
Away from the symmetry points, this formula yields back the expression Grp correspond- 
ing to complex zeros of real random polynomials: for instance, near the real axis but for 
q far from |z, we have 

G^iQ, P) ~ G'rp(P) = _^-Jy,2 (95) 

(the function Grp is shown in fig. 7 left, dashed curve). Away from the symmetry axes 
(i.e. at a distance ^ we recover a uniform density p"'^^'^(a;) ~ 1. We get the 

density on the whole torus using the same periodicity properties as in the complex even 
case. 

A separate treatment has to be made for zeros on the symmetry axes [26, 29, 30]. It 
yields a singular density on these axes: 

pT^x) = -i={i(p) + i(p-l/2)}Gr(^«) + (««p), (96) 



with Gf'(O) S ^ta.hW^) + -g^)"' (97) 



(see fig. 7 left, dotted curve). Among the zeros of iIj{z)^ the proportion situated on 
the symmetry axes is asymptotically 4 / x/tt/V, which corresponds exactly to the difference 
1 - JT2pT^^''{x)dx. 
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The description of the pair correlation function is too involved to be presented here. 
Far from the symmetry points, it reproduces the results of [28] , where it was shown that 
the function P2°^'''^(xi, ^2) takes the isotropic form (72) when xi and X2 are far from 
the symmetry axes, whereas the repulsion takes a different shape near the axes. As for 
the zeros on the symmetry axes, they repel each other as in the case of real random 
polynomials [31, eq. (5.35)], as long as they stay away from Xgy^- 

We did not compute explicitly the form factor for this ensemble, but we noticed 
qualitative modifications for the averages of the Fourier coefficients, due to the strong 
anisotropy of pi, and its singular part. Writing (pk) as a function of «; = k/y/N, the 



N 

dominant contribution is not isotropic; in case both kg and kp are even, we get the formula: 

Pk}^ ^ h + ^^^{h,FRp{Kp) + 6kpFRp{Kg)} + —Fp^intiK). (98) 

where Frp is the Fourier transform of Grp((5) + S{Q); -Fpomt is unknown analytically, it 
gives the corrections due to the symmetry points. The term in 1/ y/N is clearly anisotropic. 
We could not obtain a closed formula for Frp, but only the following hmits (see fig. 8, 
solid curve): 

Frp(«) = O(k^) as K^O (99) 
Frp{k) ~ 1 - as K ^ 00. (100) 



Notice that the dominant term oi (^(^p^J ^ ^ is of the same order 1/A^ as the form factor 
derived for the complex model (92): although we did not compute the form factor in 
the real model, we nonetheless expect the anisotropy to be visible in numerical data, 
especially for k > 1. More precisely, it seems reasonable to conjecture that the variance 

A 9 /i |9\ real / / \ rcaK 2 

Apf. — \ \Pk\ ) ~ {.y^/N J ^^^•^ different from its value in the complex model, in 

particular it should be almost isotropic. This assumption provides the following tentative 
formula for the form factor: 



real / / \ real\ 2 



In the next section, we will use this approximate formula (see fig. 6 right, dashed curve) 
to probe the anisotropy of the averaged Fourier coefficients of eigenstates for the cat map 
S and the baker's map. The fixed-A; asymptotics for \pk\ given by (101) is still oc N~^/'^: 

^\{Z)k' + ^{5^^kl + 5,X)\^ (103) 

using the asymptotics Frp(k) ~ /3ac^ as k ^ 0. (The corresponding equation for is 
easy to derive). 
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6 Results and questions about the Fourier coefficients 



We now return in greater detail to the Fourier coefficients of chaotic eigenconstellations. 
The Fourier transform of a density of N points on the torus was defined in eq. (21) as 

1 ^ 

Pfc = ^Ee"'"''"^ (104) 

Pq = 1 entirely corresponds to the normalization, but in all other respects it must be 
understood that the restriction k ^ expressly applies. The Fourier coefficients pk have 
already appeared several times above: 

1) the cquidistribution of zeros in the phase space was most simply expressed as the 
property that each pk individually tends to zero as ^ +oo; we further argued that the 
Schnirelman property could amount to a stronger statement about their rate of decrease, 
as being o{N~^) (our conjecture 1). 

2) The simple toy constellation of Sec. 4. 2. 4 carried a modulation essentially produced 
by the lowest Fourier coefficients (|/c| = 1); this largest possible wavelength was never- 
theless able to excite an extremely localized 'artificial scar'; this illustrates the idea that 
semiclassical features of Husimi densities are likely to be strongly coupled with collective 
(rather than individual) degrees of freedom of the zeros, of which Fourier coefficients are 
an example. 

3) Statistical ensembles of zeros make predictions about the moments of the Fourier co- 
efficients; especially, the model without symmetry predicts the quadratic average <^|pfcP^^ 
to be an isotropic function, universal up to the axis scales, explicitly describable for small 
or large k; the large-iV behavior (iPfcP)^ ~ ^\/C(3) \k\'^N~^^'^ followed. 

4) Eq.(12), relating the logarithm of the Husimi function to the density of zeros p{x), 
was made diagonal in the Fourier space by eq.(23), as 

hk = — ^ (except ■.ho = 0) where h{x) = e^"''' '' ■ (105) 



We will now try to build a (still partial) Fourier picture of chaotic eigenconstellations 
incorporating the previous remarks. We have just listed several arguments supporting this 
Fourier approach, but difficulties also arise which can however be worded as interesting 
problems. 

a) lack of physical cogency: any (semi)classical meaning of the Fourier space and of 
the action of classical maps upon it are presently quite unclear — save for linear (cat) 
maps — , and so are the quantum and semiclassical dynamics in this Fourier picture; 

b) besides eq. (105), another piece of the link from p{x) to the Husimi density is the 
formula H{x) — ce^^^^^ (with c = l/GM[i7]), and now this relation is local in the phase 
space, and also nonlinear; its Fourier space translation is then quite intractable, meaning 
that we cannot easily transform the Fourier coefficients of logH (the primary objects of 
interest, in our approach) into those of the Husimi density itself (which would also relate 
simply to the Fourier coefficients of the Wigner function) ; 
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c) a basic well known difficulty with Fourier transformation is that it obliterates the 
very special nature (positivity and finite parametrization) of this density p. The zeros 
being noteworthy for their 1-1 correspondence with pure states, it is specially desirable 
to fight the redundancy now reinstated by going to a countable set of Fourier coefficients, 
and this point is discussed next. 

6.1 Essential Fourier coefficients 

Problems: can one identify a minimal finite set (of size 0{N)) of Fourier coefficients pk 
that will optimally encode, and robustly restore, the coordinates of the zeros in any 
configuration? Such Fourier coefficients (dubbed 'essential'), which are in finite number, 
will then determine all the others, in a way which can also be asked. 

In one dimension this is a classic problem: the Fourier transformation (104) on a 1-d 
torus, setting Zj = e"^'^^^^, becomes {zj} {p^ = Y.jLi Zj}^ the celebrated mapping from 
sets of N points to their symmetric functions. This mapping is known to be algebraically 
invertible using the first N symmetric functions as data; the solution relies on certain 
cumulant expansions, which also happen to play a role in the semiclassical analysis of 
determinants of quantum eigenvalues. It is thus interesting to encounter a 2-d analog of 
this problem in connection with the corresponding analysis for eigenfunctions. 

This problem seems much harder than its 1-d progenitor, and we do not know if it has 
been answered. This difficulty is seen from elementary geometry (area estimates). The 
issue is to preserve and recover an amount of information consisting of 2A'" independent 
real numbers (initially: the coordinates of the zeros, essentially arbitrary). To optimize 
isotropy and robustness, we will seek this information only in radially truncated Fourier 
transforms: {pk}keDKJ = the disk {|A;| < K}; the symmetry relation p^k = Pk VA; 
(expressing the reality of p{x)) then makes each complex Fourier coefficient in the set 
count as real, while their total number can be estimated by the disk area ttK^. All in 
all, lossless encoding then requires an area of at least 2A'". It is also hard to conceive how 
the lowest 2N Fourier coefficients could generically have 'hidden' dependence relations, 
besides any obvious symmetries. So we expect the Fourier data to be sufficient starting 
from the disk radius = \fN k< with k< =^ ^2/7r ^ 0.7979. 

But let us now examine the special family of all square lattices of zeros (possibly 
rotated) for any N. The Fourier transform of such a lattice is a Dirac delta distribution 
on the dual square /c-lattice which has generators of length \/iV, so that all its Fourier 
coefficients vanish identically within any disk of area < ttN , although this area is 7r/2 
times the heuristic estimate above. The same reasoning using the optimizing lattice 
geometries in the sense of Sec. 3. 3 (lattices becoming equilateral triangular in the limit 
N — > oo) further dilates this disk of indetermination up to the area 2t^N/ -s/S. Within any 
lesser radius, then, some patterns will be totally undetectable and hence undeterminable 
(e.g., almost equilateral-triangular lattices will be undistinguishable from any of their 
translated or rotated images). Intuitively we also believe this is the worst possible case 
for the present argument. We thus expect the Fourier data to be unambiguous only from 
the disk radius K> = y/N k> with k> = ^2/\/3 1.0746. 
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So, contrary to the 1-d case where the two radii are reached simultaneously at A^, 

here there is a gap of a factor ^Jti/\/2> ^ 1.3468 from an inner radius from which we 
naively expect the information about the zeros to be generically preserved, to an outer 
radius Ky needed for the actual recovery of the zeros in some cases, and hence for their 
robust reconstruction in all cases. In-bctwccn there seems to lie a blurred zone in which 
one could perhaps invert more and more robustly by adding redundancy gradually. This 
suggests that there may not exist a unique 'best' Fourier inversion algorithm feeding on 
finite Fourier data, in contrast to the 1-d case. 

It is also interesting to consider this restoration issue for random data (complex-valued 
and without any symmetry, for simplicity). The relevant object is then the form factor 
F2{k,) with K, =^ \k\/^/N shown on fig. 6 right. For k > this form factor is virtually 
undistinguishable from being Poissonian; a gradual transition sets in around the annulus 
< K < K>}, and below «;< there lies a basin of increasingly damped expectation 
values as K — > 0: this feature makes any individual fluctuations much more significant 
against the statistical noise in this region. But precisely, the essential Fourier coefficients 
all lie within this basin, hence they are the ones statistically enhanced (especially far 
inside). By contrast, the outer Fourier coefficients, which look completely like noise 
(Poissonian) and also dominate in average size, are devoid of primary meaning since they 
ought to be (comphcated) functions of the lower ones. Moreover, the Husimi density 
itself is controlled, via h{x), by the coefficients hk of eq. (105) which have the form factor 
n'"^ F2{k) / k'^ of eq.(79), also shown on fig. 6 right. This form factor now enhances the 
coefficients in the essential region numerically, meaning that a resummation truncated 
within this region must give a fairly accurate picture of the function h{x) on a scale where 
individual zeros cannot be separated (as will be validated numerically below). 

6.2 Fourier coefficients for fixed N 
6.2.1 Averaged behavior 

We have computed the behavior as a function of k of Fourier coefficients pk quadratically 
averaged over the whole basis of eigenf unctions at fixed N, for a few classically chaotic 
models. (All of these have the parity symmetry, which makes the pk purely real.) 

We first consider a dynamical model without anti-unitary symmetry: the cat map 5", 
because the corresponding statistical form factor is a well controllable isotropic function in 
Fourier space. (This actually holds for the statistical form factor without parity symmetry, 
the only one we can fully compute, but parity-induced deviations must be small for large 
N and are indeed not visible on our plots.) Fig.9a shows the mean values of jp^p and 
along selected directions in Fourier space for this map quantized with = 149, 
restricting to even-parity eigenfunctions for definiteness. We did not compute a theoretical 
dispersion curve (involving 4-point correlation functions) , but we note that the dispersions 
stay practically uniformly within 25% of the theoretical mean value curve, however low 
the latter becomes. We conclude that quahtatively speaking, the distribution of averages 
very closely follows a random-like behavior, in a perfectly isotropic fashion. 
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The results of the same calculations are now shown (fig. 9b) for the cat map S and 
the baker's map, which both possess antiunitary symmetries. Random theory predicts 
a violation of isotropy already for (pfe)^ at large k along the directions dual to the 
symmetry lines, i.e., the two axes for the cat map vs the bisecting diagonals for the 
baker's map. Indeed, the averages (iPfcP)^^ ^-^e selectively larger in the specified singular 
directions. Isotropy improves as k 0, but the behavior stays systematically above the 
curve without symmetry. As was the case with the invariants, dispersions are noticeably 
larger for the baker than for the cat map. Still, the gross features remain qualitatively 
random-like as before. 

6.2.2 Individual behavior 

Firstly, anticipating Sec.6.3 (fig. 13), we note that within the set of individual eigenstates 
of a fixed-A^ quantum map, each given Fourier coefficient viewed in isolation will have 
large relative fiuctuations, like the invariants we studied before (fig. 3). 

We now consider the collective behavior of Fourier coefficients within individual eigen- 
constellations, by looking at fully 2-d Fourier plots of some characteristic states, earher 
portrayed by Husimi plots (figs.l, 4). Although pk and hk are interchangeable as Fourier 
representations, their plots stress quite disjoint features of the solutions. 

The striking fact about the p^-plots (fig. 10 top) is that, notwithstanding the large 
individual coefficient fluctuations, these plots neatly display a central region where most 
Fourier coefficients are heavily depressed, very much like the basin of the random model. 
This region appears for every single chaotic eigenstate we have examined in this man- 
ner, and its radius remains comparable to ^/N. The presence of this basin reflects, in 
the Fourier space, both the repulsion between the zeros and their equidistribution. We 
stress that this is now a random-like feature generically embodied in individual chaotic 
eigenconstellations (still on empirical grounds), and not in integrable ones. 

As regards the /i^-plots (flg.lO bottom), we argue that they can reveal localized peaks 
of the Husimi densities such as scars. A measure for the presence of a point scar (taken at 
the origin, up to a trivial translation) can be given by an averaged value of h{q,p) about 
that point: a relatively large and positive value will signal a scar. On the other hand, a 
convenient average over the correct width 0{N~^/^) is supplied by the partial Fourier sum 

h-K J2\k\<Khk (truncated in the disk Dk) precisely when K ^/N•. i.e., this average 
is basically controlled by the essential Fourier coefficients. We can sec indeed (fig. 11) 
that as a function of the cutoff radius K, gradually drifts towards its limiting value 
= h{0, 0) while K grows within the basin, and basically stabilizes further out (unless a zero 
happens to accidentally lie very near the origin: then h{0, 0) is very large and negative, 
but at the same time it is not a good estimator for the desired average). Moreover, the 
correct order of magnitude is often approached well inside the basin, i.e. the lower Fourier 
coefficients dominate the sum (as predicted by the shape of the form factor (79) in the 
random model). In conclusion, a scar or high density peak, taken at (or as) the origin, 
is betrayed by essential Fourier coefficients hk which, on a global trend, fluctuate away 
from their averages (hk) in the positive direction, especially for low k. (The argument 
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holds for when the coefficients are complex.) 

Fig. 12 illustrates, now over the whole phase space, three stages of the radially trun- 
cated Fourier resummation for the scarred = 128 baker cigcndcnsity of fig. lb (center), 
using the respective cutoff values K = 2, K = 9 K^), X = 12 (~ A'>). The main scar 
already begins to emerge for K = 2 0.18 V^) which, upon desymmctrization, leaves 
only 4 (real) nonzero Fourier coefficients! At the other end we may consider the Husimi 
density as recovered (in the sense of measures), reasonably well for K fa and quite 
well for K — 12 ^ K^. By contrast, the individual zeros and the associated factorized 
structure of the Husimi density are washed away by the truncations, since the Fourier 
picture is dual (complementary) to the stellar representation. The next challenge is then 
to find some effective procedure to achieve the eigenvector reconstruction itself from these 
data, i.e., to unravel the zeros which will generate this function h{x) already recovered 
on some coarse scale 0{1/VN). The shape of the form factor and the example of the 
lattice states suggest that Fourier coefficients in the transition annulus {K^ ^ 1^1 ^ 
may become critical for this purpose, all the more so when the zeros' repulsion is stronger 
and their pattern more rigid, but this example is special in that its duality is explicitly 
implementable (by the Poisson summation formula); for more disordered cases, no similar 
analytical handle is available and the restoration of the individual zeros remains an open 
problem. 

One conclusion of this discussion is that Fourier coefficients of constellations can pro- 
vide an imperfect but still effective way of reducing the number of parameters involved 
in the description of chaotic eigendensities: the number of essential Fourier coefficients is 
roughly comparable to the number of degrees of freedom of the zeros, but in practice their 
contribution to the Husimi density often decreases gradually with k to become negligible 
when or before the basin boundary is reached. 



6.3 Fixed Fourier coefficients for variable N 

We may now return to the issue of semiclassical ergodicity as expressed by the decay 
rates of fixed-A; coefficients, say = — Pfc/(7r|A;p), as oo. We present numerical 

data concerning the lower such coefficients hi^o and hi^i in graphic displays similar to 
the earlier plots of invariants (fig. 3). Fig. 13a shows the values of those coefficients for 
even-parity eigenconstellations of the cat map S, in both log-hnear and log-log plots, 
while fig. 13b shows hi^Q only for the baker's map. (Other low Fourier coefficients, not 
shown, behave likewise, whereas high {k ^ \/N) Fourier coefficients confirm the behavior 
oc A^~^/^ predicted by eqs. (101) and fig. 8.) 

The general comments made about the invariants in Sec. 3.4 carry over to these 
scalar quantities. Here, however, the conjectures linked with semiclassical ergodicity can 
be probed more sensitively through the decay rates of \hk\ for each fixed k and N ^ oo: 

one random model predicted (l^fcP)^ ^ ~ \/C(3)-^~^^^, whereas we conjectured the 
Schnirelman property to hold as long as hk = o{N~^). We can now fit various subsets of 
the above data more or less reliably with behaviors of the form cN^ for \hk\ and estimate 
the resulting values of /? (the slopes on the log-log plots). 
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Cat map S: we find that P ranges between —1.4 and —1.5 for the (quadratically) 
averaged coefficients; as for the uppermost coefficients at each N, they decrease with 
slopes reasonably steeper than —1.1. 

Baker's map: /? ^ —1.1 (resp. < —1.33) for the (quadratically) averaged coefficients 
of states whose parity is even (resp. odd, not shown); the uppermost coefficients at each 
N give overall slopes ~ —1. (resp. <. —1.2) for even- (resp. odd-) parity states: the even 
case seems to be on the borderhne of conceivable exceptions to the Schnirelman property. 

The plots also highlight the values corresponding to the eigenstates maximizing H{0, 0) 
for each (i.e., those with the strongest scars at the origin). Then the left-hand-side plots 
confirm the globally upward deviation of each Fourier coefficient of such a state from the 
linear average. The iV-dependences are erratic except, remarkably, for the baker's map 
restricted to powers of 2 where they become extremely regular, giving slopes ~ —1.4. The 
Fourier coefficients for this, or other suitable, subset(s) of N are dynamical quantities 
which could then perhaps be described asymptotically. 

7 Conclusion 

The results of this article are mainly geometrical and statistical. Indeed, the only role 
played by the ergodic nature of classical dynamics was to assert, through Schnirelman's 
theorem, the weak-* convergence (in the semiclassical limit) of the Husimi eigendensities 
towards the Liouville measure. Although such a measure-theoretic property seems a 
priori rather weak, combining it with the analytical properties of the Husimi densities 
provided direct information on both the Bargmann eigenfunctions (linearly related to the 
wavefunctions) , and their eigenconstellations (which parametrize the eigenstates optimally 
in phase space) . These analytical properties actually yielded results about the phase and 
the zeros of the Bargmann function, which seem totally immaterial quantities in the 
framework of Husimi measure theory. In the study of eigenstates, the relevant quantity 
seems to be the (scaled) logarithm h{x) = \ogH{x) of the Husimi function, linearly 
related to the density of zeros. For an integrable system, the strategy of WKB theory 
is also to build semiclassical approximations for the logarithms of eigenfunctions (see eq. 
(1)). In the chaotic case, we ended up with a WKB-like 'symbolic' description for the 
Bargmann eigenfunctions (section 4.3): 

'X^) «e^r^"'<^^'. " (106) 

Unlike its integrable counterpart (1), this formula is meaningless as a computational tool 
(the integral in the exponent is ill-defined), but it represents concisely two characteristic 
properties of chaotic eigenstates. 

First, as we pointed out in section 4.3, the above formal z-dependence of the analytic 
function \ogi/j{z) embodies the necessary asymptotic denseness of its singularities, or 
Bargmann zeros (the dislocation points of the phase of ipiz)). This phase then forms a 
fan-like pattern (see fig. 5b) which seems quite universal, in the sense that one cannot 
easily decipher the particular features of an individual eigenstate (eg. a scar, or in the 
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opposite a lattice eigenconstellation) by just looking at its phase pattern. Similarly, a 
peak in the Husimi function cannot be detected through the local density of zeros, but 
rather in giobai parameters of the whole density (Fourier coefficients, for instance). At the 
local level, zeros and amplitudes of the Husimi function seem almost totally uncorrelated. 

The second idea carried by the above formula is indeed universality. Unlike the WKB 
formula (1), where the classical action in the exponent takes the dynamics explicitly 
into account, equation (106) does not depend on the particular dynamical system we are 
studying, as long as it is fully ergodic. 

The idea of universality (partly) underlies the models of random vectors introduced 
in section 5, to which the numerical data of chaotic eigenstates are compared. As we 
recalled, such models are eigenstate counterparts of the random matrix models often used 
to mimic eigenspectra of classically chaotic systems. Until recently, they generated a 
statistical description of eigenstates in configuration space. We rather developed their 
analysis in the Bargmann-Husimi-stellar representations on the torus phase space, and 
found a generally good agreement with data from chaotic eigenstates, the comparison 
being made at different levels. 

We first studied the Husimi functions themselves, using norms and related functionals 
meant to probe the uniformity of these functions in finer topologies than the original 
weak-* involved in Schnirelman's theorem. When averaged over all eigenstates, these 
functionals agree very well with their statistical values. Individual fluctuations away from 
these average values demonstrate the existence of very ergodic or, on the contrary, more 
localized (e.g. scarred) eigenstates. 

We then explored the properties of eigenconstellations. As was already noted in 
[23, 26], the 1-point and 2-point correlation functions of the constellations fit perfectly 
with the random models, when averaged over all eigenstates. This corresponds to a uni- 
versal local interaction between the Husimi zeros, which (to this extent) resemble particles 
of a gas. We also obtained statistical predictions for the Fourier coefficients pk of the con- 
stellations (and, by linearity, for the coefficients hk of h{x)), which contain relevant global 
information. The overall shape of the form factor (i.e. a basin of small Fourier coefficients 
near the origin, surrounded by a white-noise-like sea for > v^iV) is clearly reflected at 
the level of individual eigenstates. We get a neater agreement with the statistical curves 
when averaging over all eigenstates (of a given parity). On the other hand, the main 
features of the Husimi density of a given eigenstate can be obtained by truncating the 
Fourier series of its logarithm around ^ \/N ; we therefore believe that the presence of 
a scar on an individual state, or more generally, dynamical information, is linked to the 
deviations of some low-A; (dubbed essential) Fourier coefficients from the statistical aver- 
ages. In other words, to go beyond the unprecise universal equation (106), we can hope 
to 'see a scar' directly on a finite set of coefficients hk- So far, this observation is not sup- 
ported by any dynamical explanation for the occurrence of a scar above a periodic orbit 
on a particular eigenstate (some heuristic arguments for this topic are given in [43]). The 
trouble here is that the translation of the dynamics onto the Fourier coefficients of h{x) is 
totally unclear. Otherwise, this analysis resembles the semiclassical theory for the spectral 
form factor K{t) [37]: this (now 1-d) form factor also exhibits high Fourier coefficients 
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displaying dominantly statistical universal behavior, vs low Fourier coefficients concen- 
trating the specific dynamical information (imprints of short periodic orbits through trace 
formulae). 

In spite of the aforementioned lack of a direct link between classical dynamics and 
the observed quantum phenomena, we believe some techniques and concepts presented 
here might prove fruitful to the further study of chaotic cigenstates. As explained above, 
the main new tool we used was the multiplicative properties of the Bargmann function 
induced by its holomorphy; this directly leads to the stellar representation as well as to 
the occurrence of logarithmic quantities {h{x) and the geometric mean GM[i/]). Further 
progress was made by Fourier-transforming these logarithmic quantities: for instance, 
their connections to semiclassical features proved more robust in Fourier space. The 
formula which epitomizes best these methods may be eq. (32). Similar tools were already 
used to study, on the one hand eigenfunctions of integrable systems (cf. the remark 
above), on the other hand spectral determinants, whose logarithms can be expanded (a la 
Fourier) in terms of classical periodic orbits. In these cases, the important features (the 
large values of the Husimi function, the zeros of the Bargmann eigenfunctions or spectral 
determinants) live on 1-dimensional curves. Here, we used the same tools to describe 
quantities (functions, constellations) living genuinely in 2 dimensions. The problem of 
recovery of the full constellation from a finite number of Fourier coefficients shows that 
the dimensional jump is not trivial, except maybe in very non-generic cases (for instance, 
the separable constellations considered in section 4.2.4, or lattice constellations). 

At all levels, the agreement between statistical predictions and exact eigenfunction 
data was much closer for the cat maps than for the baker's map, for which large deviations 
were repeatedly encountered. In particular, the baker's data did not exclude sequences of 
exceptional eigenstates, for which Schnirelman's property would not hold (by contrast, it 
is proven to hold for all cat eigenstates we considered). We believe that such discrepancies 
might be due to the discontinuity of the classical map, whose role had been noticed (see 
for instance [41, 40]) when studying spectral properties of the quantum map. On the other 
hand, cat maps have arithmetical properties, which make their spectra very non-generic. 
It would be interesting to find out if the pseudo-random nature of their eigenstates persists 
if we perturb the map continuously to make it 'generic'. 
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Appendices 

A Optimization of invariants for lattice states 

In this appendix we derive explicitly the geometric mean and the norm for the Husimi 
functions of lattice states, defined in sections 3.3 and 4.4. The scale invariance of these 
functionals means that we just need to derive them for the basic functions defined 
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in equation (17). We will then prove that both these functionals take absolute extremal 
values for the triangular lattice r = e''^/^ (we are indebted to M. Bauer for this proof). 

By using eq.(17) and standard formulas on theta functions [1],[2, vol. 2, chap. 13], 
one obtains the following formulae: 



GM[i7^(,)] = exp (^-^J^ \ogH^^r){q,p)dqd^ = p%T)\r,{T)\^ (107) 
Il^x(r)ll2 = ^ j^^ H^(r){q,pfdqdp = E^exp (^-tt^^ J . (108) 

In the formula above, ri{T) is Dedekind's modular form [52, p. 129] ?7(r) = e"^"^/^^ nn>i(l" 
g2i7rnr^ Curiously, these quantities are closely related to the determinants of the Lapla- 
cians A,- on T^. Indeed, we have [50, 51] 



det 

tr(e*^-) 

so Il^x(r)ll2 = tr(e*^-)|i=Q(^)/4^. 

The quantities (107,108) are obviously modular invariant, so we only need to study 
them for r in a fundamental domain of PSL(2,z) (see fig. 14). We will show in the 
following that the only extrema for both invariants on this domain are situated at the 
symmetry points r = i (B on the figure) and r = e^'^/^ (A on figure 14). 



$>(r)^|,7(r)r = GM[//,(.)]^$>(r)/2 (109) 



A.l Optimization of the geometric mean 

The equation for an extremal point of GM[if^(T-)] reads ^(r) = 4^7)) which imphes that 
3?(^(t)) = 0. If one denotes e^''^'^ = re^^, this equation reads 

nr"sin(ng) _ 

^ h _ r-^io^ne 12 ^' \^^^) 
n>l 1^ ' ^ I 

so it holds trivially for ^ = mod vr. For small r, the first term of the sum dominates 
the remainder, unless sin(^) is very small. We actually show that for r smaller than 
To ~ 0.134, the only solutions of (111) are the two points {sin(^) = 0}. Indeed, if we 
isolate the first term in (111), we obtain the inequality: 



rlsmiyjl 1.5-^ nr^smlni 



1 — r e 



ie 12 



2 



-lE ,'^' ^^^^'-T/. < E , ^ .,. |sm(g)|. (112) 

I ^ \\ — fn QinO 2 — ^ M _ Qind 2 I ^ ^ I ^ ' 



n>2 1-^ ' I n>2 



If sin(^) 7^ 0, a few simplifications yield the inequality r(l— r)^ < 4r^ — 3r^+r^, which does 
not hold for r G [0, Tq], or equivalently for 9'(r) > to ~ 0.32 (on the figure, this corresponds 
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to the dashed horizontal hne). Thus, the only extremal points with Q'(t) > tg must be on 
one of the axes {3?t = 0},{^t = 1/2}. 

It remains to study the variations of GM[H-^(^t)] along these two axes in the fundamental 
domain. On the imaginary axis, we have 

^ log GM[Hu] = -- + -^ + 27rV , 
dt ^ ^ ^ 12 At ;^^l-e-2^"*' 

which is a strictly decreasing function for t e r;^. On the other hand, this derivative 

vanishes for i = 1, since the modular transformation J — q ^ maps {t e]oo, 1]} 

onto {t e]0, 1]}, t = 1 is an extremal point; it is therefore the only one along this axis. 

To deal with the second axis {Q'(t) = 1/2}, we use the following relation, also due to 
the transformation J: 

When t describes r — 1/2^ is on the semi-circle (a, -B", ^4', /?). Fortunately, we know 
the sign of 3?(^(r)) in the domain ^{r) > to, since this quantity only vanishes along the 
symmetry axes, and in the limit Q^(t) ^ 00 we have 3f?(^(r)) ~ 27r e"^'^^^'^) sin(27r3fJ(T)). 
Therefore, ^{^{t)) takes the sign of sin(27r3fJ(r)) as long as Q^(r) > to, as indicated on the 
figure. This provides the variations of GM[iJ^(T-)] along the line {r — 1/2 + it} between 
the points a' = J (a) and /3' = J{l3): it increases from (3' to A, then decreases down to 
B', where it takes the same value as at the conjugate point B; the variations from B' to 
a' are symmetric. A direct estimation of the series logGM[i?i/2+it] shows that its sign 
does not change above 

We thus conclude that the invariant GM[if^(T-)] has only two extremal points in the 
fundamental domain of PSL(2, Z): 

-it has a global maximum at r = e"^/"^, corresponding to the equilateral triangu- 
lar lattice, i.e. the closest packing of points on the plane. There, it takes the value 
GM[i7,w3] = ^r(l/3)3. 

-it has a saddle-point at r = i, i.e. the square lattice. There, we have GM[ifi] = 

-it vanishes in the limit ^{t) — > 00. 
A. 2 Optimization of the L2 norm 

The variations of the second invariant | |-ffx(r)l I2 studied similarly. Indeed, the real 
part of the equation = reads: 

^ jA;e--^W(^'+^') sin(27rjA;3?(r)) = 0, (113) 

j,k>0 
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which yields sin(27r3?(r)) = as long as Q'(t) > ti ~ 0.3067, so the extremal points 
in the fundamental domain must also be situated on the two symmetry axes {3?(t) = 
0}, = 1/2}. The variations along the line 3?(r) = 1/2 arc obtained as above, i.e 

through a modular transformation to the half-circle {a, B', A', j3): wc obtain variations of 
opposite signs as in the former case (on this axis, | |if;^(T-)| I2 has a minimum at the point 
^4). The analysis on the axis 3?(r) = can be performed directly on the series. Precisely, 
from the formula 

and the fact that the function x 1— > e~^(l — Ax) is strictly increasing for all x > 5/4, we 
deduce that the whole derivative increases for t>l. Besides, this derivative vanishes by 
symmetry at the point t — 1. Therefore, ||i7it||2 has a single extremum on the imaginary 
axis, at t = 1. 

This second invariant is thus quite similar to the geometric mean: 

- it has an absolute minimum for the equilateral lattice ||ifei'^/3||2 = %vi^2 ■ 

- it has a saddle-point for the square lattice = H-f^xlli = ^-2^% ■ 

- it diverges as 9(t) — > 00. 



B Invariant lattices of cat maps 

In this appendix we derive a few properties of the invariant lattices of cat maps; these 
properties are used for the estimations made in section 4.4. A classical 'cat' map is given 

by a hyperbolic matrix "S* = ^ ^ d ) ^^(2, Z) (i.e. ad — be = 1, |a + ci| > 2). 

When studying the quantized map Us on TIn in Bargmann's representation [42], one 
is led to search for sublattices A of invariant under S. More precisely, A must be 
a free principal submodulc of Z%, i.e generated by a unique integer vector V: A — 
{aV mod N, a = 0, . . . , — 1}, and these N points must be different. 

A can also be considered as a lattice in Z^: it then admits the basis ^ ^ 

main task here is to estimate the minimal basis [Vi, V2] for A, i.e. the basis composed 
of the two successive shortest (for the euclidean norm) non coUinear vectors: < |Vi| < 
IV2I < ... (one can show that [Vi, ¥2] is indeed a basis for A [53, p. 83]). 

A general theorem ensures that the shortest non- vanishing vector Vi has its norm 
|Vi| < \/2N [53, p. 137]. In the following, we will prove that for invariant lattices, this 
norm also admits a lower bound of the type C\fN . This property will then imply the 
same estimates for the second basis vector V2, and a control over the size of the Dirichlet- 
Voronoi cell (in the present context, this cell can be defined as the set of all points in 
closer to the origin than to any other point of A) . 

We will study such invariant lattices in the case where they admit a generator V of 
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the form T/ = ^ j (this is possible if b and N are coprime) . The resulting sublattice A 
modulo N is invariant though 5" iff 

bk'^ + {a-d)k-c = mod N. (115) 

To estimate the minimal vector Vi, we search for the shortest vector congruent to aV 
modulo N, for a = 1, . . . ,N (we can restrict a to values in 0, . . . , N/2 by parity). Its 
square norm is |V"„p = a"^ + {(3 = min|Q;A;| mod N^. 

In all our calculations occurs an important quantity related to S, its discriminant 
D = (^)^ — 1, whose square-root is irrational because |a + (i| > 2. Equation (115) leads 
to the following constraint between a and (3: 

{b(3 - a'^-^f = a^D mod N. 

Since D is not the square of a rational, this equation cannot hold in Z. For a < ^ N/ 
it then yields the inequality 

P > l(N/MTiV - a\^\) f{a). 

We then study the variations of /(a), which depend upon the sign of the product be. 

If be > 0, f admits a minimum at ai = \^^\\f^^, where it takes a value /(cti) = 
Ci\fN . Therefore, the square norms \Va\'^ > A^min(Ci, 1/-D). 

If be < 0, f{a) decreases monotonically from a — 0, and vanishes at the point a2 — 
^jN/\be\. However, we can still bound the norms as \Va\ > C\/N , with C a positive 

constant: if a < min(Q;2/2, ^Jn/D), then jKl^ ^ fi^T > /(«2/2)^; otherwise, |KP > 
mm{ay4:,N/D). 

Prom the bounds of the shortest vector Cy/N < \Vi\ < \/2N, let us now study V2. On 
the one hand, it defines a parallelogram with Vi, of area iV, that we choose to orient as 
V1AV2 = —N (this is a choice between ±V2). On the other hand, since V2+nVi is also in A, 
we must have \Vi.V2\ < l/2|Kp. Both constraints yield \Vi\'^ < |1^2p < N/2{^-^ + ^). 

This inequality actually improves the upper bound on the minimal vector: < A'"y^4/3; 
besides, the lower bound on implies < iV(C^/4 + 1/C^). 

We can now estimate the radius of the Dirichlet-Voronoi cell centered at the origin, 
which is a hexagon except in the case V1.V2 = (a rectangle). The radius is then the 
distance of the farthest vertex from the origin. We find the uniform bound 

< :^max(5/2,C' + l/C2). 

The invariant lattice A is mapped to the constellation of an eigenstate ■(/; in TYjv through 
a rescaling and a complex conjugation; the constellation is the sublattice of Tc generated 
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by [^1,^2], with Vi = Vi/{N^/2). The above estimates then yield l^il ~ 0{\/N), and the 
modulus T = ^2/^1 is bounded in a compact set: 

y374 < 9t = A^/lyiP < and - 1/2 < SRr = < 1/2. 

r 1 1 

C Sup-norm estimates for the statistical model 

We present the derivation of the typical sup-norm of the Husimi density of a random 
complex vector, which yields the estimate listed in table 1, 3rd line. Such a quantity can- 
not be deduced from the same methods as the other invariants, but requires a completely 
different analysis. 

The proof proceeds geometrically. Instead of considering an ensemble of Gaussian 
random states (61), we constrain the statistical states to be on the unit (complex) sphere 

, with the probability density given by the canonical volume form. In Euclidean 
coordinates, the volume of this sphere reads 

Vol(5^) = / cPaoa-\ao\Y~''^o\(S'^-^) (116) 
Jd(o,i) 

= 27rnV d\{l-\aAY-^-'' 
M, Jd(oa) 



j=0 



where -D(0, 1) is the unit disk, and qq is the overlap between the normalized state iV'ia}) 
and the first vector of the basis, \qo)N,o', this overlap is related to the Fubini-Study (or 
Hermitian) distance rfps between the representatives of these two vectors in the complex 
projective plane CP^, which wc note [ip{a}] and [go]: 

|ao| = \{'^{a}\qo)\ =cosciFs(['0{a}],[go])- (117) 
Prom this, we obtain the normalized volume of the ball B^dqal^d) of radius d in CP^: 
Vol{BNi[qold)) _ Vol{{iP{ay e S"" s.t \ao\>cosd}) _ 2 



Vol(CP^) Vol(cS 



^{l-cos'df-\ (118) 



In our statistical framework, this quantity is the probability for a random normalized 
state \ip{a}) to have an overlap |(^{a}|0)| > cose? with a fixed (normalized) state 

Now, equation (64) shows that the square- norms of the torus coherent states \z)n,o 
(see eq. (3)) are given by jsf e^^^^^ in the classical limit. Therefore, in this limit the 
Husimi density H^{z, z) can be written as 

H^{z,z) ^ N\{iP\z)r,,,ra? = N cos^ d^sMM). (119) 

where [z\ is the representative of the normalized torus coherent state |-2)norm- Therefore, 
the maximum of H^{z,z) is given by the distance between [ip] and the set of coherent 
states in the projective plane CP^: {[z], z G Tc}. To estimate this distance, we first 
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select M points Zi on Tc, which are well-distributed: for instance, we consider the square 
sublattice of of side 1/-\/M, with M ~ N'^f^, for a certain power /3 > 1/2 to be selected 
later. Thus, any z ^Tq will be close to one of the zf 

^z e Tc, 3zi s.t. \z-z,\< N-^ =^ dpsiiz], [zi]) < N^/^'f^. (120) 

For a given distance d, we can estimate the proportion of states s.t. maxj | | > 
cosd, i.e. 3i, dps{[z], [zi]) < d: 

YoljUfi, B^{[z,],d)) ^ ^ Vol{B^{[z,],d)) 

Vol(CP^) - ^ Vol(CP^) ^ ^ 

< M(l-cos'd)^-^ (122) 

Using this inequality, we give a lower bound on the distance d = dj^ s.t. the proportion of 
such states in CP^ is equal to a fixed number p G (0, 1) (we use the scaling law M = N'^^): 

Vol(CP^) ^ 



, /2/31ogA^+ |logp| \ dcf .T^.x 
dN > arccos | ^/ ^^-—^ = f{p,N) (124) 

Now, if a state [ip] is not in Ui^i BN{[zi\,dN), i-e. if its distance from any Zi is greater 
than d^, the triangular inequality implies 

V^eTc, c?Fs([^],W) > dpsm,[z,])-dps{[z],[zi]) (125) 

> f{p,N) - V2^N^/^-'^ (126) 
^H^{z,z) < Ncos\f{p,N)-V27rN^/^-^) (127) 

If we fix p e (0, 1) and take P — 1, this uniform bound on behaves as 2 log A/" in the 
limit N oo. This bound is valid for a proportion (1 — p) of states in CP^^. Prom there 
we easily deduce the inequality stated in table 1 for the average sup-norm. 
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Figure captions 

Fig. la. Husimi density plots and stellar representations (plus their Voronoi tessella- 
tions) for a sample of even-parity eigenfunctions of the S-ca.t map (eq.(13)) quantized at 
N — 107, a splitting prime. The Voronoi partition of a constellation gives an idea of its 
local density of zeros (as the inverse area of each cell). 

Top 2 rows, from left to right: the most uniformly spread eigendensity; one with 
typical (random-like) values of invariants; a lattice state (the only non-real eigenstate; 
together with the conjugate lattice state they span the doubly degenerate eigenspace of 
this operator). Bottom 2 rows: the most strongly scarred eigenstate (at the fixed point 
(0,0)); another state scarred around (0,0) and along the stable and unstable hues; another 
localized eigendensity, featuring the highest H-function and lowest geometric mean. 

General caption for Husimi plots: Husimi densities are shown on a linear gray scale, 
common to all plots of each figure to preserve contrast differences (with zero as white); 
a small white triangle locates a maximum, and the white circle around it materializes a 
Planck area (equaling h — 1/N); main norms and invariants (Sec. 3. 2) are listed at top; 
zeros arc located by black stars. 

Fig. lb. Same as fig. la for a sample of strongly scarred odd-parity eigcnstatcs of 
the quantum baker's map. Leftmost plots: two members (A^ = 48, 128) of the sequence 
of those eigendensities with the highest scars above the period-2 classical orbit (through 
(1/3, 2/3)); right: an eigenstate strongly scarred along the stable and unstable hues of 
the singular fixed point (0,0). 

Fig. 2a. Some comparison Husimi densities (non localized) on the square torus; their 
constellations of zeros are now shown superimposed. Left: the singlet (A^ = 1) density 
Hy^ (cf. eq.(9)), also giving the elementary density cell of any square lattice state (cf. 
fig. 4, left); middle: an almost optimally equidistributed lattice density, for = 56 (zeros 
form near-equilateral triangles, of horizontal side 1/7 and height 1/8); right: a computer- 
produced sample of the random ensemble (without any symmetry: Sec. 5.1) for the same 
value of = 56. 

Fig. 2b. More comparison Husimi densities (locahzed) on the square torus. Left: an 
A^ = 16 plane wave; middle: an A^ = 16 coherent state (cf. eq.(3)); right: a degenerate 
density having a single zero of order A^ = 16 (white star), corresponding to the 16-th 
power of the singlet density (fig.2a, left). 

Fig. 3a Some invariant fluctuation measures of the even-parity eigenstates for the 
quantization of the S-cai map of eq.(13), with the following splitting prime values of N: 
13, 37, 59, 61, 83, 107, 109, 131, 157, 179, 227, 347, 397. The following functionals of 
the Husimi eigendensities are shown: (a) sup-norm; (b) L^-norm; (c) H-function (nega- 
tive entropy); (d) Geometric Mean (flipped upside-down). The '+' mark all individual 
values and the thick curve connects their suitable averages at fixed N (up to 131). The 
bullets pinpoint the values corresponding to the most scarred state above the fixed point 
(0,0) for each N (meaning for us the state giving the highest Husimi value i?(0, 0)); the 
dashed curve connects these values, only for visual emphasis; all curves drawn here arc 
devoid of significance outside of the computed points they interpolate between. Data 
for some comparison states of Table 1 are added: the equilateral triangle lattice value 
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(long-dashed line), the random-ensemble value (dot-dashed lines, including the typical 
deviations around the averages — except for the sup-norm, where the line is only the up- 
per bound of Table 1), and the two integrable-state examples (solid curves). 

Fig. 3b. Some invariant fluctuation measures of the eigenstates for the quantum 
baker's map, with selected even values of N up to 512. The following functionals of 
the Husimi eigendensities are shown: (a) sup-norm for even-parity states; (b) same for 
odd-parity states; (c) H-function (negative entropy) for even-parity states; (d) same for 
odd-parity states. The mark all individual values and the thick curve connects their 
suitable averages at fixed N (up to 128). The bullets mark the values corresponding to 
the most scarred state above a specific periodic point for each A^, namely: the fixed point 
(0,0) in the even-parity case, and the period-2 point (1/3, 2/3) in the odd-parity case (cf. 
the 2 leftmost plots of fig. lb); curves moreover connect two special A'"-subsequences from 
the latter set, i.e. N (dashed), and A/" = 3 x 2^= (dotted — for odd states only). The 
same comments apply to curves as in the preceding figure. 

Fig. 4. Striking effect of a long-wavelength deformation upon a constellation of zeros. 
Left: N = Husimi density of square lattice, with the periods {1/8, 1/8} (each cell 
reproduces in miniature the singlet density of fig. 2a, left). Middle: Husimi density 
of a smoothly deformed constellation (the mesh is slightly squeezed towards the edges 
and stretched towards the center to reach a sides' ratio of 0.9, almost unnoticeable on 
the constellation). Rigfit: same deformation for = 256 (meaning now a four times 
smaller Planck's constant). Notice the large variation of the sup-norms (and of the other 
functionals too); due to this, the three contrast scales are here exceptionally selected 
independently of one another. 

Fig. 5a. Illustration of the semiclassical behavior of the logarithmic derivatives 77(2;) = 
(-2) of Bargmann eigenfunctions, using color density plots to represent complex 
functions (a complex value is encoded by intensity for modulus — with zero as white — 
and hue for phase). Central plot: the classical conjugate momentum ^ as a function of z 
(from which our color encoding of complex phases can be inferred). Left plot: 7r{z) for 
a classically integrable example, an eigenstate of a quantum Hamiltonian (15-th state of 
Harper's operator — (cos27rp + cos27rg) for = 31). The lines enclose a vicinity of the 
classical energy curve, here defined to be the region where the Husimi density exceeds 
1/10-th of its maximum value; the approximation 7r{z) ~ ^ is good only nearby. Right 
plot: Tr{z) for a classically chaotic example, an eigenstate of the quantum baker's map (the 
strongly scarred eigenstate of fig. lb right, with A^ = 128); the approximation 7r{z) ~ z 
now holds throughout, save very close to Bargmann zeros, or poles of Tr{z) (the 'pimples' 
in both quantum plots). 

Fig. 5b. Plots of Bargmann eigenfunctions il'{z), using intensity to encode the 
Husimi density and color to encode the phase. Left: the 15-th eigenstate of the quantum 
(N — 31) Harper Hamiltonian (same as in preceding figure), specifically with a binary 
coding of the modulus (threshold = 1/10-th of maximum Husimi density) letting the 
phases show through everywhere in the classically forbidden regions. Middle: an (odd- 
parity) eigenstate of the quantum 5'-cat map for the same A^. Rigfit: the eigenstate of 
the quantum (A^ = 128) baker's map of preceding figure and fig. lb right. A universal 
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radial fan-like pattern for the quantum phases, having the regular wavevector distribu- 
tion {kg, kp) ~ ttN X (p, —q) prevails almost everywhere in the case of classically chaotic 
systems (barring dislocation points at the Bargmann zeros), and only about the classical 
orbit in integrable cases. 

Fig. 6. Left: the rescaled pair-correlation function g{r) between zeros of a random 
Bargmann function with no symmetry (eq. (72)). The shape (short-distance repulsion 
followed by a small bump) is similar to the corresponding correlation for particles in a 
dilute gas [32]. The scaling r = x/tt/V \5z\ indicates that zeros repel each other up to a 
distance ~ and are uncorrelated beyond. Right: various rescaled form factors. 

Solid curve: rescaled form factor F2{k) for the density of zeros p{x) in the model with 
no symmetry, eq. (76). Dot-dashed curve: rescaled form factor for the logarithm h{x) 
(eq.(22)). These form factors also hold in the semiclassical limit for the model of even- 
parity complex random states (cf. eq. (87)), and conjecturally for the model of even 
real random states as well, in all directions except the axes (eq. (101)). Dashed curve: 
conjectured rescaled form factor along the two symmetry axes for the density of zeros of 
even real random states. These curves are compared to numerical data in section 6. The 
two dotted vertical lines represent the rescaled radii «;< and k>, described in section 6.1. 

Fig. 7. Left: average densities of zeros, for various statistical models. Solid curve: 
average density of zeros for the model of even-parity complex random states, as a function 
of the rescaled distance from a symmetry point (eq. (86)). The drop of the density in 
the vicinity of the origin can be interpreted as a repulsion between Zi and —Zi, whereas 
the subsequent bump ensures the global normalization of pi. The dotted and dashed 
curves concern the model of real even random states. Dashed curve: density of complex 
zeros transversally to an axis of symmetry, far from symmetry points, so that parity 
corrections are negligible (this curve is the vertical cut {Q = 3} in fig. 7 right). Dotted 
curve: renormalized density of singular zeros (e.g. real zeros) away from a symmetry 
point (eq. (96)). Right: Contour plot of the average density of complex zeros near a 
symmetry point for the model of even-parity real random states, as a function of the 
rescaled variables {Q,P) = VttN x {q,p). The contour level spacing is 0.03. We only 
show the first quadrant {Q > 0,P > 0}, the others being obtained by reflection w.r.t. the 
axes. The density vanishes along the axes and has the bulk value 1 far from the symmetry 
axes. The maximum of Gi (?» 1.084) hes at the point Q = P ^ 1.042. As long as \Q\ or 
\P\ > 2, a good approximation is given by the formula for random polynomials (95). The 
cut of this plot along the axis {Q — 3} is given by the dotted curve of fig. 7 left. 

Fig. 8. Solid curve: rescaled plot of the dominant (oc N~^^'^) term in the expansion 
of the average Fourier coefficients along the symmetry axes for the model of even-parity 
real random states (eq. (98)). This dominant term is of the same order as the quadratic 
average (101), so it should be visible on the numerical data. Dashed curve: rescaled 
average Fourier coefficients away from a symmetry point for the model of even-parity 
complex random states (86). This average is of order N~^, whereas the quadratic average 
(92) is of order N'^/^. 

Fig. 9a. Square averages of the Fourier coefficients of the density of zeros (left), and of 
the logarithmic Husimi density (right), related by eq.(23), for the even-parity eigenstates 
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of the ^'-cat map with N = 149. The abscissa carries two scales: the modulus of the 
Fourier index (in boldface), and its rescaled value \k\/^/N. The Fourier coefficients are 
computed for k along selected directions from the origin, labeled by distinctive symbols: 
+ horizontal fc-axis, x vertical fc-axis, A first diagonal, y second diagonal. The data are 
compared to the statistical form factors in the complex random state model (see fig. 6 
right, and eq. (92)). 

Fig. 9b. Same as fig. 9a, but for the even-parity eigenstates of the 5'-cat map (top) 
and baker map (bottom) and also including the data along the direction (2, 1) (marked 
with o). Form factors of the real random state model are given both for a generic direction 
(dot-dashed curves) and along a symmetry axis (long-dashed line) (cq. 101). Due to the 
large deviations from the statistical curves in the baker's case, we use varying ordinate 
scales. 

Fig. 10. Two-dimensional density plots for the Fourier coefficients pk {0 k e 1?) 

(top row) and hk = —pk/i^k^ (bottom row) of some constellations. Positive values are 
in red, negative values in blue (zero is white). The two circles mark the inner and outer 

radii, = ^2N/'k and K-^ = \J2N/ respectively. Density scales arc independently 
selected, so we indicate the largest absolute value appearing in each plot ('Max'). 

Left: most ergodic N — 107 5'-cat eigenstate of fig. la, top left; middle: most scarred 
N — 107 S-c&t eigenstate of fig. la, bottom left; right: N — deformed lattice state of 
fig. 4, middle. 

Top row: the two leftmost p^-plots (Max ^ 0.357) exemplify two generic features of 
chaotic eigenstates: the strong squeeze of the coefficient values in the central basin, and 
their noisy look outside; the rightmost p^-plot (Max ^ 0.956) is an example where the 
essential coefficients begin to rise only in the transition region. 

Bottom row: the Fourier coefficients of h{x) = N~^logH{x) show distributions of 
widths ~ V^. The leftmost /i^-plot (Max ^ 0.00134) corresponds to a Husimi density 
with a very moderate value H{0) fa 1.397) and shows a good balance between positive 
(red) and negative (blue) values of h^. The middle (Max Ki 0.00161) and right (Max 

0.00837) /ifc-plots correspond to densities very peaked at x — 0, and show a marked 
overall dominance of the positive (red) values of hk- 

Fig. 11. Partial sums x J2\k\<K as functions of (integer) K for the six N = 107 
5'-cat map eigenstates of fig. la (each curve is labeled by the value of the corresponding 
eigenangle in units of 27r marked at the hmiting value logi^(O) on the right edge, except 
for one state having \ogH{0) —9.3). The inner and outer radii and are marked 
by the dotted vertical lines. The limiting behaviors are seen to set in gradually before 
K reaches (the oscillations are a Gibbs phenomenon, which would diminish under 
smoother truncation schemes). 

Fig. 12. Coarse-grained Husimi density of the = 128 baker eigenstate of fig. lb 
(middle), obtained through partial Fourier summations of its logarithm, truncated at 
radius K = 2 (left), 9 ~ (middle), 12 fsi (right). (This is not a smoothing of the 
Husimi density itself by convolution.) The resulting densities are normalized a posteriori 
(by brute force), then the other approximate invariant functionals are computed. 

Fig. 13a. The low Fourier coefficients hi^ (top row) and hi^i (bottom row) of the even- 



61 



parity eigenstates for the quantization of the ^-cat map. In left column: log-hnear scale 
plots of N^l'^hk, in the same general conventions as fig. 3a. In right column: log-log plots 
of the solid curve is the quadratic average at fixed A^, against the dot-dashed straight 
line of the random model (without symmetry) which shows the slope —3/2. As in fig. 3a, 
the values corresponding to the most scarred state above the fixed point (0,0) for each N 
are marked by bullets connected by the dashed curve. 

Fig. 13b. Same as fig. 13a, but for the even-parity eigenstates of the quantum baker's 
map, and only for the lowest Fourier coefficient hi^. As in fig. 3b (left column), the bullets 
mark the values corresponding to the most scarred state above the hxed point (0,0) for 
each A^, and the dashed curve connects these values above the special subsequence A^ = 2*^. 

Fig. 14. Poincare's upper-half-plane {Q'(t) > 0}. The fundamental domain of 
PSL(2, Z) is bounded by the 3 pieces of hyperbolic geodesies drawn as thick curves. The 
horizontal dashed line corresponds to 9(t) — to'- above this line, ^{r]'/r]) has the sign of 
sin(27r3?(r)), as marked. The modular transformation J maps the semi-circle (a, 5" , A', j3) 
onto the vertical axis (a', -B', A, /?'). Large arrows indicate directions of increasing values 
of GM[i/^(^)] along the semi-circles and the vertical axes 3?(r) = 0, ±1/2. The maxima 
of GM[Hy^{r)\ are the points A, A! and their modular images. 
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sup H GM sup H GM sup H GM 

12.2484 1.5015 0.4457 0.5797 7.5477 1.4780 0.4865 0.4898 6.2266 1.5116 0.5457 0.4569 
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sup H GM sup H GM sup H GM 

3.9959 1.2673 0.2971 0.6491 6.1711 1.3913 0.4058 0.5545 2.4056 1.3043 0.3871 0.5289 




arXiv:chao-dyn/9711016 vl 15 Nov 1997 




arXiv:chao-dyn/9711016 vl 15 Nov 1997 




Baker P = + Baker P = - 
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